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ABSTRACT 

We present a Chandra, Suzaku and Rosat study of the hot Intra Group Medium (IGrM) of the 
relaxed fossil group/ poor cluster RXJ 1159+5531. This group exhibits an advantageous combination 
of flat surface brightness profile, high luminosity and optimal distance, allowing the gas to be detected 
out to the virial radius (R v ; r = i?ios =1100 kpc) in a single Suzaku pointing, while the complementary 
Chandra data reveal a round morphology and relaxed IGrM image down to kpc scales. We measure the 
IGrM entropy profile over ~3 orders of magnitude in radius, including 3 data bins beyond ~ O.5R200 
that have good azimuthal coverage (>30%). We find no evidence that the profile flattens at large 
scales (>R5oo), and when corrected for the enclosed gas fraction, the entropy profile is very close to 
the predictions from self-similar structure formation simulations, as seen in massive clusters. Within 
Rvir, we measure a baryon fraction of 0.17 ± 0.02, consistent with the Cosmological value. These 
results are in sharp contrast to the gas behaviour at large scales recently reported in the Virgo and 
Perseus clusters, and indicate that substantial gas clumping cannot be ubiquitous near R V i r , at least 
in highly evolved (fossil) groups. 

Subject headings: Cosmology: dark matter — Xrays: galaxies: clusters — Galaxies: groups: individual: 
RXJ1159+5531— Galaxies: ISM 



1. INTRODUCTION 

Galaxy groups (defined here as bound systems with 
virial masses, M V i r , in the range ~ 10 13 -10 14 Mq) are 
essential ingredients in the assembly of structure within 
the Universe. Locally ~30% of galaxies ar e found in 
group s that are less massive than 10 14 M Q (Eke et al. 
2004). Systems of this mass range have been implicated 



an entropy profile that rises with radius, r, such that 
S oc JjiZ—j where Tn is a characteristi c tempe rature 
dTozzi fc Normanl [200l| ; |Voit ct al.| ^005| ; |Kaisei| |l986| ) 



Observations of galaxy clusters in fact reveal lower-mass 
objects to be systematically offset upward s from this re- 
lation, but approaching it by ^Rsoo (e.g. Ponman et al 



2003 ; Pratt et al. 2010 ). This indicates nongravitational 



as key sites in b oth efficient star formation (Springel & 
Hernquist 200£ ), and the morphological transformation 
ot galaxies ( Zabludoff & Mulchacy 1998] ). Typically the 
dominant baryonic component is, however, an extended, 
hot gas halo, that c an contain > 70% o f the baryons at 
M 5 oo^5 x 10 13 M Q dGiodini ct al] |2009| fl Such a mas- 



energy injection, likely a consequence of preheating of 
the gas in filaments prior to accret ion, or feedback due to 
star formation or a central AGN ( Voit fc Ponman 2003 ; 
& Donahue 200E; Ponman et 



Voit 



2003) 



balance between these different processes s 



The fine 
rould affect 



sive hot gas halo can make them both luminous X-ray 
sources, and, potentially, targets for detection in large 
numbers by future Sunyaev-Zeldovich surveys (Haiman 
|st al.||200l| ). 

X-ray observations provide the best means for study- 
ing emission from the hot gas in groups. Of particular 



inter est is the entropy proxy, 5* 



-2/3 



kT (where n e is 



both the shape and normalization of the entropy proxy 
profile, making it a potentially powerful diag nostic (e.g. 
Borgani et al. 2005; McCarthy et al. [2010 ), especially 
at the group scale, where these effects should be com- 
paratively more important. Recent Chandra and XMM 
observations of relaxed groups have revealed complex en- 
tropy profiles that flatten at lar ge (^O.lR.^nn) and smal l 
scales, albeit with some scatter ( Gastaldello et al] 2007a • 
Mahdavi et al.||2005t j Finogucnov ct al. 2007 1: Humphrey 



and " is the temperature), which is related to tire spe 
cific entropy through a logarithm and constant ottset. it 
purely gravitational processes shape the energetics of the 
gas, self-similar structure formation simulations predict 
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the e lectron number density, k is Boltzmann's constant ct al] 2008 "Sun et a 



son ct al. 2009; Flohic et al 2011), possibly indicating 



that both star formation and AGN feedback arc im por 
tant ( |Voit fc Donahue] |2005| ; |McCarthy et al.|[2010| ). 

Feedback also shapes the enclosed gas fraction (f gas ) 
profiles of g roups and c luster s, which are fou nd to rise 
with radiusflAllen et al] |2002| ; |Vikhlinm ct~aT] [2006| ; [GOT] ; 
sun ct al. 2009), but with groups containing a system- 
atically s maller fraction of their g as at small scale s (e.g. 
Fig 11 of Humphrey et al. 2011, hereafter |H11 ). As- 
suming they can be converted into a reliable estimate 
for the total (initial) baryon fraction in the cluster (ft,), 
f gas measureme nts are a powerful cosmological tool, u sed 



either directly flWhitc et al] |1993| ; [Allen ct al] [2002P , or 
in employing the gas mass as a virial mass proxy (e.g. 



2 



Voevodkin & Vikhlinin 2004). Massive clusters are pre- 
ferred in this analysis since i gas measured at these scales 
should be closer to fb than in poor clusters or g roups. 
Still, system s with masses as low as ~ 2 x 10 14 Mfn flAllen 
et al. [2008[ ), or even lower (I Voevodkin fc Vikhlinin] |2004| ) 



atic errors in the George et al. background treatment. In 
the outermost b i ns, th e dep rojection procedu r e ado pted 



are routinely used. Trans 



ating into fb the 



by Urban et al 
pends sensitive 



values 



measured at scales far smaller than the virial radius, 
Rvir (typically at ~R25oo), involves a number of assump- 
tions that have yet to be robustly ve rified, especially in 
lower mass halos (e.g. Arnaud 2005). Intri guinglv, for 
most of the g roup-scal e obi ects studied by Gastaldcllo 



( [2011 ) and aimionescu et al. (2011) de- 
y on the correct modelling of projected 
emission fro m regions outside the field of view. Further- 
more, while [Simioncscu et al attempted to mitigate the 
scattered light contamination from the cluster core, it is 
unclear how sensitive their results were to this correction. 



were 



et al| ( 2007b , hereafter G07 ), extrapolating f gas outside 



the field of view yielded global fb constra ints consistent 



with the Universal baryon fraction (0.17: Dunkley et al, 
[2009| ; |Komatsu ct"aT[|201l|) , in accord with the idea that 



X-ray Tlrigh t groups are baryonically closed ( Mathews 



ay 

et al[| |2005| ). Nevertheless, this extrapolation was sub 
ject to significant systematic uncertainty. 

To date, most observations of the gas in groups and 



Similarly, the XMM-Newton profiles of Urban et al 
sensitive to the treatment of the background. 

In this paper, we present a joint Chandra and Suzaku 
study (carefully cross-checked with the archival Rosat 
data) of the very relaxed galaxy group RXJ 1159+5531, 
allowing, for the first time in a system with a mass as 
low as ~ 10 14 Mq, the gas to be traced to scales as 
large as the virial radius, i?ios- A single, giant ellipti 
cal galaxy dominates t he stellar light, makin g it a proto 
typical "fossil group" (Vikhlinin et al. 1999] ) and i mply 



clustsrs have been restricted to within ^Rsoo, and typ- 
ically much smaller scales are attained. For example, 
R500 was only reached for 11 of the 43 systems stu died 



in the current largest group sample (Sun et al. [2009 ). In 
groups at these scales, f gas is still only ~0.07, in contrast 
to ~0.11 for clusters. Given the low, stable background 
for the Suzaku XIS instrument, recent work has begun to 
push measurements of the ICM in clusters out to ~R200i 
or beyond, but a coherent picture has not yet emerged. 
While some studies have found c o nsistency with mode l 
predictions (Rciprich et al. 2009; Hoshino et al._2010). 



devia tions from hydrostatic equilibrium (Bautz et al 
[2009), temperat ure asymmetries associat ed with large- 



scale structure ( Kawaharada et al. 201C ) and, in three 
systems, an unexpected flattening al the entropy proxy 
profile ou tside ~R Rn o (PKS0745-1 91: |George et alj2009 



Perse us: Bimionescu et al 
|201l[) , have bee n seen 
et all! (see also Nagai 



2011; Virgo: 



ing a highly evol ved (i.e. relaxed) dynamical state (Pon 
man ct aT| 1994). The optimal combination of distance 
(z = 0.081), mass (R v i r = 12') and high surface bright- 
ness (suc h that gas is already detect ed to ~Rfinn wit h 
Chandra: |Vikhlinin et~aT[ |2006| ; |G07| ; |Sun et al]|200g| ), 
make it possible to measure the gas to R V h- in a single, 
modest (85 ks) Suzaku pointing. 

The group was observed with the X-ray centroid 
slightly offset (by 5') from the Suzaku optical axis, to en- 
able coverage out to ~12', while minimizing the compli- 
cations of stray-light from having the X-ray peak outside 
the field of view. Although this configuration limited the 
number of radial bins we can study, by combining the 
Chandra and Suzaku data, we were able to achieve ~3 
spatial bins outside ~ O.5R200J with >30% azimuthal 
coverage, in comparison with ~6 bins, and ~15% cover- 
age in the ~3 times lo nger Suzaku exposure of Perseus 
(Simionescu et al. 2011). In conjunction with the archival 



Chandra data, we were able to measure the gas properties 
over almost three orders of magnitude in radius. 
Previous studies of 



Chan dra data (Vikhlinin et al 



2006 




G07 


; Bun et al. 



2009 ), have reported conflicting parametcrizations of the 



rban ct al 
jTmioncscu ct al. | and |Urlmi i 
Lau 2011) attributed the en 
tropy flattening, and associated over-estimate of f gas , to 
putative dumpiness of the ICM in the outskirts of the 
cluster, leading to a systematically biased gas density 
measurement. 

Given the lack of a consistent story in the outskirts 
of clusters, the ubiquity of a clumpy ICM remains to 
be determi ned. The azimuthal tem perature variations 
in Perseus ( Simioncscu ct al. 2011 ) and the large-scale 
asym metrie s in the X-ray image of Virgo (Bohringei 
et al| |l994 ) indicate that these systems are not relaxed 
at large scales, consistent with ongoing formation. This 
could give rise to deviations from sphericity (complicat- 
ing the deprojection) or local hydrostatic equilibrium 
(hence an underestimate of the gravitating mass, and 

an over-estimate of f gas ), and local distortions in the en- outlined in Hll, which enable finer control of s ystematic 



gravitating mass pro file. Adopting the po pular Navarro- 
Frenk- White ( NFW : [Navarro et al] |l997|) profile (plus, 
in the case of |G07 , a baryonic component), M500 in- 
ferred from these studies has varied significantly from 
~ 6 x 10 13 -10 14 M Q , and the corresponding NFW con- 
centration pa ram eter from ~1.7-5.6, which remains a 
puzzle (see § |4.2| ). The addition of new data at large 
scales should help pin down the mass profile more pre- 
cisely, pa rticularly if the scale radius were as high as 



found by Vikhlinin et al 
employed the 



(400 kpc). In this paper, we 
forward-fitting" mass analysis techniques 



trop y profile. On account of the proximity of these clus- uncertainties tha n more traditional methods (Buote fc 



ters, only a small fraction of the outer annuli were im- 
aged niT4ielic - sTu7lT^^ 

effects would be seen in azimuthally averaged profiles. 
An additional concern is the need to minimize sources 
of systemati c uncertainty in this background-dominate d 
regime (e.g., IRciprich ct al.|[2009| , |Bautz et al. | |2009| , [mil ). 



Humphrey 2011a) 



We assumed a flat cosmology with Ho — 70km s _1 
and tt\ = 0.7. We adopted i?ios as the virial radius 
(Rvir) , based on the approximation of Bryan & Norman 
(|1998D for the redshift of RXJ 1159+5531. Unless other- 



Eckert et al.|~(2011) demonstrated that the Rosat PSPC 



surface brightness profile of PKS0745-191 disag rees at 



7.7-cr with t he Suzaku density profile inferred by George 



et al. ( |2009| ). They attributed the discrepancy to system- 



wise stated, all error-bars represent 1-a confidence limits 
(which, for our Bayesian analysis, implies the marginal- 
ized region of parameter space within which the inte- 
grated probability is 68%). 
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2. DATA REDUCTION AND ANALYSIS 
2.1. Chandra 

The region of sky containing RXJ 1159+5531 was im- 
aged by the ACIS instrument aboard Chandra on two 
separate occasions. We consider here only the deep data 
taken in the ACIS-S configuration (Observation ID 4964; 
beginning on Feb 11 2004). A shallower ACIST observa- 
tion was also available, but to simplify the analysis (in 
particular, the background modelling), we chose not to 
include it in our stu dy. T he data-reduction was carried 
out as described in Hll, using the CIA 04.1 and Hea- 
soft 6.8 software suites, in conjunction with the Chandra 
calibration database (Caldb) version 4.1.2. Briefly, the 
data were reprocessed from the "level 1" events files, fol- 
lowing the standard data reduction threads!]. Periods of 
high background were identified by eye in the lightcurve 
from a low surface-brightness region of the CCDs and 
data from these intervals were excised, leaving a total 
exposure of 75 ks. Point sources were detected in the 0.3- 
7.0 keV image with the wavdetect CIAO task, which was 
supplied a 1.7 keV exposure map to minimize spurious 
detections at chip boundaries. The detection threshold 
(10~ 6 ) guaranteed < 1 spurious source per CCD. All de- 
tected sources were confirmed visually, and appropriate 
elliptical regions containing ~ 99% of the source photons 
were generated. 

In Fig H we show a smoothed, flat-fielded Chandra im- 
age, having removed th e poin t sources with the algorithm 
outlined in Fang et al. ( 200E ). The image was smoothed 
with a Gaussian kernel, the width of which varied with 
distance from the nominal X-ray ccntroid according to 
an arbitrary powerlaw, ranging from ~1" at the centre 
of the image to ~1' at its edge. The image is smooth 
and very round, consistent with the relaxed morphology 
expected for a fossil system. To search for more subtle 
structure, we used dedicated software to fit an elliptical 
beta model (with constant ellipticity) to the central 2'- 
wide portion of the unsmoothed (flat-fielded) image^. In 
Fig[l], we plot (data— model) 2 /model, corresponding (ap- 
proximately) to the x 2 residuals from this fit. To bring 
out the structure, we smoothed this image with a Gaus- 
sian kernel of width 3 pixels. There is weak evidence of 
a small (~4"), coherent structure in the residuals within 
the central (^10 kpc) region, which may imply a small 
depression in the surface brightness (by ~20-30%). The 
formal significance of this feature depends sensitively on 
prior information (in particular, the region of the im- 
age oyer wh ich one searches for structures; e.g. Kaastra 
tal|p006| ). Still, even if the feature is real, it should 
not give rise to a significant error i n the recovered, az- 
imuth a llv averaged mass prof ile (see Buote & Humphrey 



E 



at Churazov et al 



2008| ) . Aside from this modest 
feature, the overall lack of significant, coherent residuals 
indicates that the X-ray image is very relaxed. 

Spectra were extracted in a series of contiguous, con- 
centric annuli centred at the X-ray centroid. The widths 
of the annuli were chosen to contain approximately the 
same number of background-subtracted counts, while en- 
suring sufficient photons for useful spectral analysis. The 

7 http:/ /cxc. harvard.edu/ciao/threads/index. html 

8 We obtained a best-fitting major axis core radius of 2", =0.51 
and axis ratio of 0.9. 



resulting annuli had widths larger than ^3", which is 
sufficient to prevent spectral mixing between adjacent 
annuli on account of the finite spatial resolution of the 
mirrors. Data in the vicinity of point sources and chip 
gaps were excluded. We extracted spectra from all the 
active chips (excluding S4, which suffers from noise). Ap- 
propriate count-weighted spectral response matrices were 
generated for each annulus with the standard CIA O tools 
mkwarf and mkacisrmf . Representative spectra, without 
background subtraction, are shown in Fig ^. 

Spectral-fitting was carried out in the energy band 
0.5-7.0 keV, using Xspec vers. 12. 5. In, by minimizing 
the C-statistic to mitigate biases that arise (even in the 
high-count regime) when using the sta ndard y 2 approx- 



imatio ns for Poisson-distributed data ( Humphrey et al 
[2009b| ) . To aid convergence, we rebirmed the spectra 
to ensure at least 20 photons per bin. The data in all 
annuli were fitted simultaneously, to enable the source 
and background com ponen ts t o be modelled at the sam e 
time. In keeping with Hll and pastaldello et al. ( 2007t| ), 
we modelled the projected (rather than the deprojected) 
source emission in each annulus as coming from a single 
APEC plasma model with variable a bundances, modified 
by fo reground Galactic absorption (Dickey & Lockman 
1990 ) | We allowed the total abundance of Fe (Zp c ) and 
the abundance ratios with respect to Fe of the elements 
O, Ne, Mg, Si and Ni to vary in each annulus. The 
abundance of He, and the abu ndance ratios of the o ther 
elements were fixed to 1 Solar ( Asplund et al. 2004 ). To 
improve S/N, we tied Zp e between the outermost two 
annuli, and constrained the abundance ratios to be the 
same in all annuli. 

To account for emission from undetected LMXBs in 
the central galaxy, we included an additional 7.3 keV 
bremsstrahlung component. Since the number of X-ray 
point sou rces is approximately prop ortional to the stellar 
light (e.g. Humphrey fe Buotq 200Sj ) , the relative normal- 
ization of this component between each annulus was fixed 
to match the relative K-band luminosity in the match- 
ing regions, which we measured from the 2MASS image. 
This component is only important in the very central 
region (<20 kpc) of the system. 

To accommodate the background, we included ad- 
ditional spectral components in our fits. Specifi- 
cally, we included two (unabsorbed) APEC components 
(kT=0.07 keV and 0.2 keV ) and an (absorbed) power - 



2004) 



law component (r = 1.41; |De Luca fe Molendi 
The normalization of each component within eaci 
lus was assumed to scale with the extraction area, but 
the total normalizations were fitted freely. We discuss 
the possible impact of an additional, " Solar wind charge 
exchange" background component in § |6.2. 1 . To account 



for the instrumental background, we included a number 
of Gaussian lines and a broken powerlaw model, which 
were not folded through the ARF. We included sepa- 
rate instrumental components for the front- and back- 
illuminated chips, and assumed that the normalization 
of each component scaled with the area of the extraction 
annulus which overlapped the appropriate chips. The 
normalization of each component, and the shape of the 
instrumental components, were allowed to fit freely. We 



We discuss the results from a deprojected analysis in 





Right ascension 

Fig. 1. — Suzaku XISO (top right) and Chandra (bottom left) images of RXJ 1159+5531 . The apparent asymmetries in the XISO image 
are due to the asymmetric point spread function of the telescope, rather than the intrinsic shape of the source. The Chandra image has 
been cleaned of point sources and mildly smoothed. The smoothing scale varied from ~1" at the smallest scales, to ~1.1' (100 kpc) at the 
outer part of the image (see text). At the top left, we show a "residual significance" image (see text) of the centre of the group, indicating 
deviations from a smooth model fit to the (Chandra) X-ray isophotes. There is no obvious large-scale feature in this map, indicating the 
system is largely relaxed. 

included two Gaussian lines (at 1.77 and ~2.2 keV), the 
intrinsic widths of which were fixed to zero. The ener- 
gies of the ~2.2 keV lines were allowed to fit freely, as 
were as the normalizations of all the components. To 
verify the fit had not become trapped in a local min- 
imum, we explored the local parameter space by step- 
ping individual parameters over a range centred around 
the best-fitting value (analogous to computing error-bars 
with the algorithm of Cash 1976). The covariance matrix 
(which contains the error bars) was com puted via the effi 



density to be constrained (this is similar to the outermost 
annuli of NGC 5044, studied by |Buote et alj ^004|) . 



cient M onte Carlo technique outlined in Humphrey et al 
( p006| ), and we carried out 250 error simulations. 

The best-fitting models are shown in Fig ^] for a rep- 
resentative selection of spectra. While the instrumen- 
tal background is clearly significant at high energies 
(> 2 keV), only in the outermost annulus does the source 
signal fall below the background level. Nevertheless, 
given the optimal temperature of the gas (~1 keV), the 
Fe L-shell peak is still visible as a small "bump" in the 
spectrum at ~0.9 keV, enabling the gas temperature and 



2.2. Suzaku 

The region of sky containing RXJ 1159+5531 was im- 
aged by Suzaku beginning on Feb 5 2009 (observation ID 
804051010), with three of the XIS units operating. Data- 
reduction was performed using the Heasoft 6.8 software 
suite, in conjunction with the XIS calibration database 
(Caldb) version 20090925. To ensure up-to-date cali- 
bration, the unscreened data were re-pipelined with the 
aepipeline task and analysed following the standard 
data-reduction guidelines^. Since the data for each in- 
strument were divided into differently telemetered events 
file formats, we converted the "5 x 5" formatted data into 
"3 x 3" format, and merged them with the "3 x 3" events 
files. The lightcurve of each instrument was examined 
for periods of anomalously high background, but no sig- 

10 http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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Fig. 2. — Representative Chandra and Suzaku XIS1 spectra for RXJ 1159+5531, shown (Chandra) without background subtraction, 
or (Suzaku) with the instrumental background subtracted. In addition to the data, we show the best-fitting model, folded through the 
instrumental response (solid black line), along with the decomposition of this model into its various components. For the Chandra data, 
we show the hot gas contribution (solid red line), the composite emission from X-ray binaries (dash-dot magenta line), the instrumental 
background (dotted orange) and the cosmic X-ray background (dashed purple line). For the XIS1 data, we show the hot gas emission from 
each annulus as solid lines (red, green, blue and orange indicating, respectively, the 0-2', 2-5', 5-9'and 9-13' apertures, respectively) , sky 
background (using the same colour scheme as for Chandra). For Chandra, the background is dominated by the instrumental component, 
whereas for Suzaku, which has a lower instrumental background but is less able to resolve the cosmic component into individual point 
sources, the cosmic component dominates. 



nificant amount of data was found to be contaminated 
in this way, leaving 85 ks of total "cleaned" exposure 
time. 85 ks total exposure time survived flare cleaning. 
In Fig 0, we show the 0.5-7.0 keV image for the XIS0 
detector, excluding data in the vicinity of the calibration 
sources. By visual inspection of the images for all three 
active detectors, we found only 1 bright point-source (co- 
inciding with one of the calibration sources in the XIS0 
image). In subsequent analysis, we excluded a circular 
region of 2.5' radius, centred on this source (which should 



eliminate > 90% of the source photons from contaminat- 
ing any spectra). 

Since RXJ 1159+5531 was slightly offset from the cen- 
tre of the field of view, it was possible to extract spectra 
out to scales of ~13' from the single, pointed observation. 
Therefore, spectra were extracted in four concentric an- 
nuli (0-2', 2-5', 5-9' and 9-13'), centred at the nominal 
position of RXJ 1159+5531 in the field of view of each 
instrument. Due to the large field of view of Suzaku, 
our spectral extraction regions actually achieve ~55-60% 
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azimuthal coverage in the 5-9' aperture, and ^27% in 
the 9-13' region. Data in the vicinity of the calibration 
sources and the identified point sources were excluded. 
For each spectrum, we generated an associated redistri- 
bution matrix file (RMF) using the xisrmf gen tool and 
an estimate of the instrumental background with the 
xisnxbgen task. Ancillary response files (ARFs) were 
generated for each s pectrum with the xissimarf gen tool 
([shisaki et al. 2007), which models the telescope's optics 
through ray-tracing. Since the point-spread function of 
the telescope is very large (~2' half power diameter), 
even with such large apertures it is necessary to account 
for spectral mixing between each annul us. W e did this 



by employing the algorithm described in Hll 



We modelled the spectrum in each annulus as the sum 
of source plus sky background components. The instru- 
mental background (whi ch is not the dominant back- 
ground component; Hll) was subtracted directly, us- 
ing the model generated with the xisnxbgen task. For 
the source, we included an APEC component, for which 
the Fe abundance and abundance ratios of other species 
were allowed to vary (as for the Chandra data). The 
abundance ratios were tied between all annuli, and Zp c 
was tied between the outermost two annuli. In the 
central bin, we also included a 7.3 keV bremsstrahlung 
component, to account for possible LMXBs in the cen- 
tral galaxy. For this component, we obtained a to- 
tal X-ray luminosity per unit K-band optical light of 
Lx.lmxb/Lk = 2.7 ± 3.9 x 10 29 er<? s" 1 ^ 1 , in good 
agreement with the Chandra result of 6.0 ± 1.9 x 
10 29 er<? s~ 1 Lq 1 . This self-consistency gives us confi- 
dence in our treatment of the background (described 
below), and the lack of bright source contamination in 
the Suzaku field. This ratio is also marginally consis- 
tent with_thc_jn£aiiscah^ determined from local galax- 
ies ( Humphrey et al. 2008 ). 

The Chandra data reveals a strong temperature gradi- 
ent within the central ^2'. Since we can generally ap- 
proximate the integrated spectrum of a region in which 
there is a tem perature gradient by the sum of two APEC 
components (Buotc 1999; Buote et al. 2003| ), we added 
a second APEC component in the central bin. Although 
we did not consider the fitted temperature or density in 
this bin in our subsequent analysis, it was still necessary 
to fit a reasonably accurate model to the data, as some 
fraction of the photons from this region are scattered 
into the surrounding annuli. The spectra from all in- 
struments were fitted simultaneously, allowing additional 
multiplicative constants in the fit to enable the relative 
normalization of the model for each instrument to vary 
with respect to the XIS 0. 

We found that this model was able to fit the data 
wellp]. We show in Fig || representative spectra for the 
XIS1 instrument, showing various model components. 
As is immediately clear, mixing between annuli is sig- 
nificant. Furthermore, the data are clearly background- 
dominated in the outermost annulus (while the source 
and background are comparable in the 6-10' region). 

11 Since the C-statistic is not easily interpretable as a goodncss- 
of-fit, we also .performed a \ 2 fit Although, in general, such fits are 



more biased (Humphrey et al 



Nevertheless, given the ^keV temperature of the gas, 
the Fe L-shell peak is still clearly visible over the cos- 
mic X-ray background component (§ ||), and so the gas 
properties can be measured with reasonable precision. 

2.3. Rosat 

RXJ 1159+5531 was observed serendipitously by Rosat 
in the outskirts (~17' off-axis; ~2.5' inside the support- 
ring structure) of a pointed PSPC observation (obser- 
vation ID 700055n00). Given its relatively poor spa- 
tial resolution and limited energy pass-band, we only re- 
quired to measure the Ros at s urface brightness profile, 



centred at the source (see § 4.5 for more details). We do 



not expect much extended emission from the principal 
target of the observation, an unrelated, foreground Syl 
galaxy (NGC3998). Therefore, we extracted the surface 
brightness profile of RXJ 1159+5531 out to scales of ^25' 
(~ 2R V i r ). To prepare the data, we followed t he standard 



Rosa t data-reduction recipe (see, for example Fang et al 
2009). A full field of view image with 15" pixel size was 



generated in the 0.42-2.0 keV band, and the correspond- 
ing exposure maps were generated with the pcexpmap 
task. Point sources were detected with the wavdetect 
CI AO task to search for structure at scales of 1, 2, 4 and 
8 pixels, and supplied with the exposure map to mini- 
mize spurious detections at the image boundaries. The 
detection threshold was set to 10 -6 , implying < 0.1 spu- 
rious detections per image, and all point sources were 
visually confirmed. We used the ao and castpart Hea- 
soft tasks to generate, respectively, the scattered Solar 
X-ray background and the particle background contribu- 
tion, which were then subtracted from the image. The 
surface brightness profile was generated in concentric an- 
nuli from the flat-fielded image, excluding data from the 
vicinity of any detected point source and, to maximize 
the ability to identify and exclude point sources, all pho- 
tons not inside the shadow of the inner PSPC support 
ring. We discuss the surface brightness analysis in detail 
in 8 fl5. 



3. SPECTRAL FITTING RESULTS 

3.1. Metal abundances 

In Fig ||, we show the projected abundance pro- 
file, which is centrally peaked, as expected for a re- 
laxed galaxy group and consistent with a picture in 
which metal enrichment is facilitated by m ass-loss from 



2009b)), the best-fitting parameters 



were very close to those obtained with the C-statistic, and the fit 
was formally very good (x 2 /dof=1787/1811). 



stars and super novae in the central galaxy ( jMathews fe 
Brighcnti 2003). There is good agreement between the 
profiles measured outside ~100 kpc with Chandra and 
Suzaku, supporting our treatment of these data. To ef- 
fect a comparison within the central 2', where there is a 
strong abundance gradient revealed by Chandra, but only 
a single Suzaku data-bin, we extracted a single Chan- 
dra spectrum for this region, and fitted it with the same 
model as the Suzaku data. The best-fitting Zp c (shown 
on Fig [s]) was consistent (within ~1.5-cr) with the Suzaku 
measurement. The good agreement over the entire radial 
range is encouraging, and suggests that our treatment of 
the spectral mixing between the Suzaku annuli is approx- 
imately correct. 

We note the slight dip in the central Chandra bin, sim- 
ilar to features that have been reported in some other 
galaxy groups and clusters. Still, if hot gas components 
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Fig. 3. — Projected Fe abundance profile, measured with Chandra 
(triangles) and Suzaku (stars). The dashed line indicates the Chan- 
dra abundance measured in same region as the central Suzaku bin. 
Note the overall good agreement between Chandra and Suzaku. 
We also overlay (crosses) the deprojected abundance profile mea- 
sured with Chandra, including a second temperature component in 
the innermost bin. Note that the central dip seen in the projected 
profile is absent in this case. 




Ne Mg 
Element 

Fig. 4. — Ahund anrp ratios with rpsp prt to Fe, expressed relative 
to the Solar ratios (Usplund et al. 20041). Chandra data are marked 
with triangles, and suzaKu data witn stars. The solid line is the 
best-fit model (x 2 /dof=5.7/8) where the enrichment entr ies from 
SNIa and_SNII, assuming the WDD2 SNIa yields frnrn 
et al.| ( |l997b| ) and the SNII yields from [Nfomoto et al 
The 



enrichment fraction is 0.82 ±0.05 



For reference, we also 



show the best fits with the W7 (green dashed line) and WDD1 (blue 
dotted line) SNIa yields, neither of which are formally acceptable. 



with a range of different temperatures are found along 
the line-of-sight to the central annulus, the measured 
Zp e may b e systematically under - estimated due to the 



"Fe bias" flBuote k Fabian| |1998| ; |Buotcj |2000b|) . Such 
a situation could arise cither due to projection effects, 
or if there is a strong temperature gradient within that 
annulus. To investi gate this, we c arried out a depro- 
jection analysis (see 



Hll 



or ij 6.4 of this paper, for a 
detailed description of the deprojection procedure), and, 
furthermore, included an additional APEC component 
within the cental bin, with abundances tied to those of 
the other component. Since the deprojection procedure 
tends to increase the error-bars on the measured quanti- 
ties, we found it necessary to fix the abundance outside 



~400 kpc to 0.2 Solar (consistent with the projected re- 
sults). We found that the second hot gas component was 
formally required in the central bin (the improvement in 
the C-statistic was 13.7 for a difference of 2 degrees of 
freedom when it was added); the two temperatures in 
this bin were 1.29 ± 0.23 keV and 0.72 ± 0.17keV, consis- 
tent with a continuing temperature decline in the group's 
centre, although it may also reflect deficiencies in the de- 
projection procedure. 

The deprojected Zp c profile is shown in Fig [|, and it 
also exhibits a centrally-peaked shape and, significantly, 
no evidence of a central abundance dip. This strongly 
suggests that the feature seen in the central bin of the 
projected data is simply an artefact of the Fe bias. 

In addition to the abundance profile, the data also pro- 
vide interesting constraints on the abundance ratios with 
respect to Fe of various species. These are summarized 
in Fig ^; we find excellent agreement between Chandra 
and Suzaku. Overlaid on this figure, we show the best- 
fitting abundance ratio patterns predicted by combin- 
ing the metal yields for type la and type II supernovae, 
through which Fe and the a elements are primarily pro 



cessed. We adopted the SNII yields from Nomoto et al 
(1997a), an d explored three differ ent theoretical SNIa 



yields from Nomoto et al. (1997a), specifically the so 
called "W7" , "WDD1" and "WDD2" models. We found 
the WDD2 model fits the data well, provided 82 ± 5% 
of the Fe was synthesized in type la supernovae. This 
is consistent with measurements in the central parts of 
other galaxy groups and clus ters (Humphrey k Buotc 
and references therein; Werner et al 



2006 



3r a 



review) . 

3.2. Temperature and density profiles 

The projected gas temperature and densityF] profiles 
are shown in Fig g| for Chandra and Suzaku. These pro- 
files span almost 3 orders of magnitude in radius. The 
outermost Suzaku d ata-bin reaches ^1200 kpc , which is 
~Rvi r , as found by Gastaldello et al. (2007b); see also 



2l There is exce 



lent agreement between the results 
for the two satellites, giving us confidence in our treat- 
ment of the mixing between Suzaku annuli, and our treat- 
ment of the background components, which differ slightly 
between both satellites. Since two APEC components 
were employed in the central Suzaku bin, we do not show 
the gas density or temperature for that bin. The Chan- 
dra temperature prof ile agrees well with that found by 
Vikhlinin et al. ( 2006| ). In Fig iL we also show the depro- 
jected temperature profile (see § \A for a description of 



how this was obtained) , which is in good agreement with 
the deprojected p r ofile re ported for the Chandra data by 



Gastaldello et al 
As discussed in 



there is evidence that the gas 
temperature may continue to fall in the central bin, pos- 
sibly requiring a second APEC component to be added 
to the spectrum in this region. Nevertheless, the results 
for the single temperature fit can still be interpreted in 
our mass-fitting analysis, provided care is taken to com- 
pute an appropriately weighted average gas density and 

12 We define the projected density in any given annulus as the 
mean gas density if all the emission measured in the annulus origi- 
nates from a region defined by the intersection of a cylindrical shell 
and a (concentric) spherical shell, both of which have the same in- 
ner and outer radii as the annulus. 
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Fig. 5. — Radial temperature (top panels) and density (bottom panels) profiles for RXJ If 59+5531. We show the projected Chandra 
profiles in the left column (triangles), the projected Suzaku profiles in the right column (stars), and the deprojected profiles in the right 
column. Overlaid are the best hydrostatic model fits to each dataset, which match the data very well and, for the deprojected profiles, 
we show in grey the f-cr confidence region for the model. We exclude the central Suzaku bin from the fitting (although it was used for 
abundance determination) since two temperature components were required, fn account for the strong temperature gradient. In the upper 



( 2006 ) from the Chandra data (light blue crosses), which. 



left panel, we show the projected temperature profile measured by Vikhlinin ct al 

agrees well with our data. In the bottom right panel, we show (light blue dashed line) "the deprojected density profile found by Vikhlinin 
et al.ll(bOOa), which is also in reasonable agreement with our data. 



temperature in that region (see Hll for more details). 
In any case, we found that our results were relatively in- 
sensitive t o th e inclusion or omission of that particular 
ann ulus Qj j6.6T). 



4, MASS MODELLING 



4.1. Method 

We translated the density and temperature profiles 
into mass constraints using the entropy-based "forward- 
f itting" techniqu e de veloped in our recent papers 
(iHumphrev et al.||2008| , [2009a| ; |mT| ; [Buote fc Humphrey 
[20114 for a review) . Briefly, given parametrized pro- 
files of "entropy" (S=kTn e 2 ^ 3 , where n e is the electron 
density) and gravitating mass (excluding the gas mass, 
which is computed self-consistently), plus the gas density 
at some canonical radius (for which we used 10 kpc), the 
three-dimensional temperature and density profiles can 
be calculated, under the hydrostatic equilibrium approxi- 
mation and fitted to the data. This model assumes spher- 
ical symmetry, which is a standard assumption in X-ray 
hydrostatic modelling, even when the X-ray isophotes 
are not perfectly round, since deviations from spheric- 
ity likely only contribute a very small error (~few per- 
cent) into the recovered ma ss profil e and barvon fraction 



tion of Prugniel fc Simicn] ( |1997 ). We fixed the effec- 
tive radius and luminosity of this component to the K- 
band values (9.8 kpc, and 1.0 x 10 12 L Q at the distance 
to RXJ 1159 +5531). inferred from 2MASS ([Gastaldello 
et al | |2007b| ), and allowed the K-band M/L ratio to be 
fitted freely. We additionally included a supermassive 
black hole, with mass fixed at 2.4 x 10 9 M©, based on the 
Lk of the galaxy and the blac k-hole mass versus V- band 
bulge luminosity relation of flGultckin ct al. 2009 ), as- 



suming Lv/Lk = 0.24, which is typical for an old stellar 
population with Solar abundances^. 

To fit the three-dimensional entropy profile we assumed 
a model comprising a broken powerlaw, plus a constant; 
such a model is reasonably successful at reproducing the 



wide radial range ( 


Humphrev et al. 


2009a 




Jetha ct al. 


2007; Finoguenov ct al 


2007| Gastaldello et al. 


2007a; 


Bun et al. 2009|). 


In order to provide more flexibility 



data-points in the fit), we also allow an additional break 
in the entropy profile at large radius. The normalization 
of the powerlaw and constant components, the radius of 
the breaks and the slopes above and below them were 
allowed to fi t freely. 



( Buote fc Humphrey 2011b; see also Piffarctti ct al. 2003 



For the gravitating mass model, we assumed an NFW 



( Navarro et al. 1997| ) dark matter halo, the virial mass 
and concentration of which were free fit parameters, plus 
a model for the stellar mass. Since the projected stel- 
lar light of the central galaxy is known to be well-fi tted 
by a de Vaucouleurs model ( Vikhlinin et al. 199E), for 
the stellar mass component, we adopted a deprojected 
de Vaucouleurs model, using the analytical approxima- 



Following Humphrey ct al. (2008), we solved the equa- 
tion of hydrostatic equilibrium to determine the gas prop- 
erties as a function of radius, from 10 pc to a large radius 
outside the field of view. For the latter, we adopted twice 

13 We note that the updated black hole mass versus K-band 
luminosity relation of |Graham| (2007, their eqn 14) gives a slightly 
smaller mass for the black hole (1.8 X 1O 9 M0). Small differences 
in the black hole mass will not affect our results, however, as the 
black hole mass is only ~2% of the total mass at the effective centre 
of the innermost Chandra bin. 
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the virial radius of the system defined by ignoring the 
baryonic components (which is slightly smaller than the 
true virial r adiu s of the system), but explore alternative 
choices in § 3.4. For any arbitrary set of mass and en- 



tropy model parameters, it is not always possible to find 
a physical solution to this equation over the full radial 
range. Such models were therefore rejected as unphysical 
during parameter space exploration. In order to compare 
to the observed projected density and temperature pro- 
files, we projected the three-dimensional gas density and 
te mperature, using a procedure similar to that described 



Gastaldello et al. 



fl2007b|f 4 |. This involves computing 
an emission-weighted projected mean temperature and 
density in each radial bin, that can be compared directly 
to the dataQ 

To compare the model to the data, we used the % 2 
statistic, fitting simultaneously the Chandra and Suzaku 
temperature and density profiles, with the central bin 
of the Suzaku data excluded from our fits. Correlations 
between the density errors were simply implemented by 
adopting a form for the x 2 s tatistic which incorporates 
the covariance matrix (e.g. Gould 2002 Parame- 
ter space was explored with a Bayesian Monte Carlo 
method. Spe cifically, we used version 2.7 of the Multi- 
Nest codePl flFeroz fc Hobson| [2008| ; [Feroz et alj ^009| ). 



Since the choice of priors is nontrivial, we followed con- 
vention in cycling through a selection of different priors 
and assessing thei r impact on our results (we discuss this 
in detail in § 6~5| ). Initially, however, we adopted flat 
priors on the logarithm of the dark matter mass (over 
the range 10 12 < M < 10 16 M Q ), the logarithm of the 
dark matter halo concentration, cdm (over the range 
1 < cdm < 100), the logarithm of the gas density at the 
canonical radius, the stellar M/Lk ratio, and the various 
entropy parameters. For a m ore detailed descr i ption o f 
the modelling procedure, s ee Humphrey et al. (2009a); 
Buote fc Humphrey! (|2011a| ). 



4.2. Mass profile 

We found the model could fit the density and tem- 
perature data well, with a best-fitting x 2 /dof=14.4/15, 
even when taking into account t he co variance between 
adjacent density data-points (se e |H11|) . The best- fitting 
models are overlaid in Fig [|. In Fig |H| we also show the 
range of possible three-dimensional temperature and den- 
sity profiles predicted by our model that are consistent 
with the projected data. In Table |l|, we tabulate the 
best-fitting values and marginalized confidence regions 

14 In computing the plasma emissivity term, we approximated 
the true three dimensional abundance profile with the projected 
abundance profile (Fig H); since the projected and deprojected 
profiles do not differ significantly, this should be sufficient for our 
eaent purposes, but we explore this question in more detail in 



6.7 



10 We note that, for systems with a strong temperature gradient 
and kT > 3 keV the emissio n- weighted temperature may be biased 
2004). Nevertheless, we find little evidence 



low (Mazzotta et al 



of any strong bias when we fit the deprojected temperature and 
density data, which should be much less sensitive to this effect. 
Thi£_suggests that the impact of this effect is not significant here 

lu By default we only consider correlations between the density 
data, but we investigate introducing a more complete covariance 
computation in § |6.9| and find that the results are not significantly 
affected. 

17 http:/ /www. mrao.cam.ac.uk/software/multinest/ 
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Fig. 6.— Radial mass profile of RXJ 1159+5531. The solid 
(black) line indicates the total enclosed mass (and the grey shaded 
region indicates the 1-cr error in the total mass distribution), the 
dashed (red) line indicates the stellar mass, the dotted (blue) line 
is the dark matter, and the dash-dot (magenta) line is the gas mass 
contribution. Overlaid are a set of data-ppints derived from, a more 
traditional analysis method described in Humphrey et al. (2009a). 
We stress that the model is not fitted to these data, but is derived 
independently from the temperature and density data; neverthe- 
less, the agreement between the different approaches: is good We 
also show the mass profile found by Vikhlinin et ah] ( 2009) from the 
Chandra data (light blue dash-dot-dot line J, wlich overall agrees 
with our results, although is higher at both large and small radii. 

for each mass parameter of interest. The best-fitting 
radial mass distribution is shown in Fig kj, along with 
the relative contributions of each of the different mass 
model components. Overlaid are a series of mass data 
points derived from fitting the deprojected data using the 
"traditional smoothed invers ion" mass modelling method 
( Buote fc Humphrey 2011a ), which is more su bjec t to 



systematic uncer t ainties (for more details see § 6.4 and 
Humphrey et al. 2009a ). The agreement is very good, 
indicating that the resulting mass profile is not overly 
sensitive to the analysis method. 

We obtain a marginalized R v i r =1100 ± 31 kpc, 
which compares to an outer Suzaku annulus spanning 
820-1190 kpc (thus having an "effective centre" at 
^1000 kpc). Therefore, we confirm that we are able 
to reach R V i r with these data. In Fig [?], we show 
the relation between the concentration of the gravitating 
mass, c v ir, and the virial mass. Mw - To be consistent 
with our past work (Buote et al. 2007), the M v ; r and R V j r 



are derived from the distribution of the total gravitating 
mass, not just the dark matter, and the concentration, 
c v ; r is defined as the ratio of R v ; r to the characteristic 
scale of the DM halo. We compare our results with con- 
centration and mass constraints for RXJ 1159+5531 re- 
ported in the literature (and based only on the Chandra 
data) at different overdensities. We fi nd excellent agree- 
ment w ith the M v i r -c v ir results from Gastaldello et al 
( 2007b| ) ; in particular the agreement with the virial mass 
and concentration obtained from that paper is interest- 
ing, given that the Gastaldello et al. ( E007b|) results at 
this overdensity were significantly extrapolated. We do, 
ho wever, find a sys tema tic offset at R^n n from the work 



of pun et alj ( 2009j ) and |Vikhlinin et al | ( 2006[ ), who em- 
ployed essentially th e same "smoothed inversio n" mass 
modelling approach (Buote & Humphrey 2011a) in each 
work, and obtained slightly higher masses and lower con- 
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TABLE 1 

Mass results and error budget 



Test 


M„/L K 
MqLq" 1 


log M 2 500 
[Mq] 


log C2500 


log M500 
[M Q ] 


log C500 


log M vir 
[Mq] 


log c vir 


Marginalized 
Best-fit 


0.54 ±0.10 
(0.51) 


13 4Q+0-03 
lo '* 3 -0 .02 
(13.49) 


0.44 ±0.07 
(0.48) 


13.77 ±0.04 
(13.75) 


0.76 ± 0.06 
(0.80) 


-,0 07+o.or, 
10 - 31 -0.04 

(13.95) 


1.05 ±0.06 
(1.09) 



ADM profile 
AAC 

ABackground 

ASWCX 

AStray light 

ARmax 
A3d 

AFit priors 
AStars 
AWeighting 
APSF 

AInstrumcnt 

ASpectral 

ADistance 

AEntropy 

ACovariance 



+0.19 v _ 011 

-0.209 (±°:i° 

+0.02 /+0.10"\ 
-0.06 ^-0.13 J 

-0.193 (±g;ii 
— 0.005 (±0.10) 

±S:S? (±°- 10 ) 

+0.05 (±0.09) 

+0.03 



£01 (±0-11) 

-0.031 (r o o ° 

0.09 J 



+0.09 



±0.01 (±0.11) 

18:81 (±o-ii) 

+0.05 (±0.10) 

±8:ia (±°- 13 ) 
+0.01 (±8 : « 

+0.02 f+omS 

-0.03 1-0.15 J 



±0.003 (±0.03) -0.015 (lo os) 



+0 - 01 (±0.03) 



0.02 
0.027 (±0.04) 

0.01 



+0.02 

v -0.03 

torn (±0-03) 
-0.049 (±0.04) 
±0.01 (±0.03) 

±0 (loTa) 
+0.03 (±0.03) 

+0.01 (ir 2 3 ) 

±0.03 (±0.03) 
±0-03 (t° : g) 

+0.01 f+0.04\ 
-0.03 ^-0.02 J 

±0.005 (±0.03) 

±S:SS (±0-03) 



±S:S? (±°- 07 ) 
±0.10 (±g : s) 

-0.006 (±0.07) 

±0M (±°- 07 ) 
±0.02 (±0.08) 
±o°ol (±0-08) 
±0.08 (±0.10) 
-0.082 (±0.07) 
-0.012 (±0.07) 

toil* (±0-0 8 ) 
-0.019 (±0.07) 

±0.11 (±0.07) 

-0.007 (±0.07) 

±0.04 f+ - ^ 



±0-06 (1 ;° 2 3 

±0.008 (±0.04) 

±S:Ss (±°- 04 ) 

-0.052 (±0.05) 

±0.003 (±0.04) 

-0.018 (+g;«) 
-0.053 (±0.05) 
±0.01 (±0.04) 

-0.020 (+g;» 
±0.04 (±0.04) 
+0.01 (±0.05) 

±S:Si (±0-05) 

±0.02 (±0.04) 

+0.03 f+0.05 
-0.08 ^-0.03 

-0.003 (±0.04) 
-0.028 (±0.04) 



-0.013 



^-0.08 J 

tao? (±°- 07 ) 
±0-09 (t° °t) 
-0.006 (±0.06) 

±0.01 (±0.06) 
±0.02 (±0.07) 

18:8s (±°- 07 ) 
±0.07 (±0.09) 

-0.074 (±0.06) 

-0.011 (±0.06) 

toicS (±0-07) 
-0.018 (±0.06) 

±0.10 (±0.07) 

-0.007 (±0.06) 

03 ( + 006 



±0.21 (±0.02) 

±0.02 (tS:£) 
18:82 (±°- 04 ) 

-0.047 (±0.05) 

±0.006 (±0.04) 

±0.01 (±0.04) 
-0.040 (±0.05) 

±S:8§ (±0.05) 
-0.022 (+8:8* 

+0.05 (±0.04) 
+0.02 (±0.05) 

±8:S (±0.05) 

+0.02 (±0.04) 

18:81" (±0.04) 

±0.006 (±0.04) 
-0.017 (±0.04) 



-0.017 (±0.07) 

±8:8? (±o.o7) 
±0.09 (18;8I) 

-0.007 (±0.06) 

±0.01 (±0.06) 
±0.02 (±0.07) 
±o.ol (±0-07) 

±0.05 (t8:5?) 

-0.069 (±0.06) 
-0.012 (±0.06) 

±8:81 (±o.o7) 

-0.016 (±0.06) 
±0.09 (±0.06) 
-0.006 (±0.06) 
±0.03 (±0.06) 



Note. — Marginalized values and 1-er confidence regions for the stellar mass-to-light (M*/Lk) ratio and the enclosed mass and concentration measured 
at various overdensities. Since the best-fitting parameters need not be identical to the marginalized values, we also list the best-fitting values for each 
parameter (in parentheses). In addition to the statistical errors, we also show estimates of th« error budget from possible sources of systematic 
uncertainty. We consider a range of different systematic effects, which are described in detail in § pt specifically wc evaluate the effect of the choice of 
dark matter halo model (ADM), adiabatic contraction (AAC), treatment of the background (ABackground) and the Solar wind charge exchange X-ray 
component (ASWCX), stray light (AStray light) maximum radius used in the projection calculation (AR max ), deprojection (A3d), priors on the model 
parameters (AFit priors), treatment of the stellar light (AStars), removing the cmissivity correction (AWeighting), the effect of errors in our treatment of 
spectral mixing due to the PSF of Suzaku (APSF), the X-ray detectors used (AInstrumcnt), spectral fitting choices (ASpectral), distance uncertainties 
(ADistance), the parameterization of the entropy model (AEntropy), and covariance between the temperature and density data-points (ACovariance). 
We list the change in the marginalized value of each parameter for every test and, in parentheses, the statistical uncertainty on the parameter determined 
from the test. Note that the systematic error estimates should not in general be added in quadrature with the statistical error. 
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Fig. 7. — Derived relation between mass and concentration for 
RXJ 1159±5531 (contours). We show the results for several over- 
densities (virial=108, in black; 500 in red; 2500 in blue). We 
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centrations than our best-fitting values. We discuss the 
possible origin of these discrepancies in § 7.2. 



4.3. Gas and baryon fraction constraints 

In Table 0, we tabulate the best-fitting values and 
marginalized confidence regions for the gas and baryon 
fractions of the system at different overdensities. In per- 
forming this calculation, we included the mass in baryons 
of the central galaxy and the hot gas (inferred from our 
fit), and also folded in two additional, uncounted reser- 
voirs of baryons, intra-clustcr light and additional galax- 
ies in the group. Vikhlinin et al. (1999) found that ~25% 
of the V-band stellar light associated with galaxies was 
in non-central galaxies. We assumed that this also holds 
in the K-band, and adopted a K-band M/L ratio of 1 for 
these galaxies. Furthermore, for the virial mass of the 
system, we expect as much as ~2 times the stellar mass 
of the central galaxy wi ll be bound up in intra-cluster 
light (Purcell et al. 2007). If wc make the reasonable as- 
sumption that these uncounted baryons are distributed 
in the same way as the dark matter (this is appr oximately 
true o f the additional galaxy light reported by Vikhlinin 
3t al.), the total mass inferred from our modelling (which 
did not explicitly include these components) will be cor- 
rect. We discuss the sensitivity of o ur results to our treat- 
ment of these components in § 6.6. In practice, however, 
these uncounted components account for only ~10% of 
the baryon budget, since most of the baryons are tied up 



11 




Radius (kpc) 

Fig. 8. — Left: Radial distribution of the enclosed gas fraction 
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TABLE 2 

Baryon fraction results and error budget 



Test 


/g,2500 


/g,500 


fg,vir 


/i),2500 


/b,500 


fb ,vir 


Marginalized 
Best-fit 


0.048 ± 0.001 
(0.048) 


0.087 ±0.005 
(0.089) 


0.14 ±0.02 
(0.149) 


n 19 +0.009 

u - iz -o.oi 
(0.117) 


124+ 007 
(0.127) 


0.17 ±0.02 
(0.174) 



ADM profile 
AAC 

ABackground 
ASWCX 
AStray light 
ARmax 
A3d 

AFit priors 
AStars 

AWeighting 
APSF 

AInstrument 

ASpectral 

ADistancc 

AEntropy 

ACovariancc 



±0.002 (±0.001) 

, f+0.002^ 
iU ^-0.001 J 

±0.001 (±0.001) 
+0.001 (±g;gS 

±0 (±0.001) 
±0 (±0.001) 
+0.004 (±0.002) 

±0 (±0.002) 

, f +0.001 
=cu ^-o.ooi 

i n /+0.002 
=cu ^-0.001 

±0 (±0.002) 

±0.002 (±r ° 2 2 ) 

-0.003 (±0.001) 

±g:Si (±o.oo2) 

±0 (±0.002) 
+n (+0.002 \ 
±U I -0.001 J 



-0.005 (±0.004) 
±0 (±0.005) 

t°;°g* (±0.006) 

±0.003 (±0.005) 
-0.002 (±0.005) 

±0.002 (ir ° 6 6 ) 

+0-01 (t°oZl) 
-0.002 (±0.005) 
±0.002 (±0.005) 

-0.002 (±0.005) 
-0.002 (±0.005) 
-0.003 (±0.01) 
-0.007 (±0.005) 
too? (±0.01) 
±0 (±0.005) 
±0.004 (+g;°°g) 



-0.044 (±0.01) 
-0.003 (±0.02) 

±S:S (±°- 02 ) 

+0.02 
-0.007 (±0.02) 
±0.004 (±0.02) 
+0.03 (±0.02) 

±g:SSi (±0.02) 

+0.009 (±0.02) 

-0.015 (±°;°2 
-0.010 (±0.02) 
-0.028 (+°;°§) 
-0.020 (±0.02) 

to.ol (±0.02) 
-0.001 (±0.02) 

±S:8Si (±0-02) 



±0.03 
-0.022 

+0.00 

-o.oo 
-0.021 
-0.003 
-0.001 
+0.02 

+0.00 
-0.00 
+0.03 
-0.04 

+0.003 
-0.005 
-0.006 
-0.004 
+0.003 
-0.003 
-0.002 



(±0.01) 
(±0.01) 
(±0.01) 

/+0.02^ 

l-o.oij 
(±0.01) 

(±0.01) 

(±0.01) 

(±0.01) 

(±0.01) 

(±0.01) 
(±0.01) 
(±0.01) 
(±0.01) 
(±0.01) 
(±0.01) 

f+0.0l\ 

K-o-oi) 



+0.006 
0.008 
+0.006 
0.007 

(±0.01) 

/+0.01 
I -0.01 



±0 
-0.012 

±0.00 
-0.005 

-0.003 £±0.01) 
-0.001 (±0.01) 
±0.02 



h0.01> 

-0.01 ; 



+0.00 
-0.00 

±0.02 

-0.003 
-0.003 
-0.005 

+0.00 
-0.01 
+0.01 
-0.00 

±0 

+0.00 
-0.00 



(±0.01) 
(±0.01) 

(±0.01) 
(±0.01) 
(±0.01) 
(±0.01) 
(±0.01) 
(±0.01) 
(±0.01) 



+0.02 

-o.oi 
±0.02 
-0.005 
±0.009 



-0.046 (±0.01) 
-0.009 (±0.02) 
(±0.02) 

+0.0l\ 
-0.02 J 

±0.02) 
(+o.oi\ 

^-0.02 J 

+0.04 (±0.02) 

+ -°ofo (±0-02) 
±0.01 (±0.02) 

-0.016 (±0.02) 
-0.006 (±0.02) 
-0.026 (±° : °i) 
-0.017 (±0.02) 
(±0.02) 
±0 (±0.02) 
+0.01 (±0.02) 



+0.02 
-0.01 



Note. — Marginalized values and l-cr confidence regions for the gas fraction (/ 9i a) and baryon fraction (fb,A) measured at vari- 
ous overdensities (A). We q4so provide the best-fitting parameters in parentheses, and a breakdown of possible sources of systematic 
uncertainty, following Table hi. We find that fb is reasonably robust to most sources of systematic uncertainty, especially within R.2500- 



in the hot gas halo. 

In Fig pi we display the radial dependence of enclosed 
bary on fraction, which rises modestly from f D ^0.1 at 



100 epe to 0.17 ± 0.02 by R V i r , in excellent a greement et al 



with t he (Josmological baryon tra ction (0.17: Dunklcy 



et al[j gOOgt |Komatsu et al.| p011| ). We stress tK ere is 
no extrapolation in this measurement of % at the virial 
radius, since the Suzaku data reach the virial radius in 
this system. This behavio ur is strikingly similar to the 
isolated Galaxy NGC 720 (Hll), which has a virial mass 
^2 orders of magnitude smaller than RXJ 1159+5531. 



In Fig ||, we also show the radial distribution of the en- 
closed gas fraction, which rises steeply with ra dius (as is 
typically ob served in groups and clusters, e.g. Vikhlinin 



|2006 ). For comparison, we overlay the predic - 
tions ol recent numerical simulations ( Young et al. 2011 ) , 
which systematically under-estimate the true gas frac- 
tion. At small scales (or higher overdensities) the low f gas 
indicates that gas has either been bound up into stars or 
"pushed out" to large radii by feedback. However, the 
approximate baryonic closure of the system suggests that 
little gas has been evacuated completely from the system 
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Fig. 9. — Contours of f gas versus mass at different over-densities. 
The magenta (dash-dot) line indicates the Universal baryon frac- 
tion (annotated "baryonic closure"). We overlay results from the 
literature. 

in such a process. That gas has been "pushed out" in this 
way is reflected in the differential gas fraction (i.e. the gas 
density divided by the total mass density at a given ra- 
dius), which actually exceeds the Cosmological baryon 
fraction outside ~500 kpc (Fig || right panel). In all, 
we find that ^65% of the gas within R v j r actually lies 
outside R500, illustrating the importance of the Suzaku 
data for directly measuring it. 

In Fig ^|, we compare our enclosed f gas constraints at 
different overdensities with values reported in the liter- 
ature. Once agai n , our m easurements agree well with 
|Gastaldello et al~ (|2007b), but the re ar e some modest 
discre pancies with Sun ct al.| @009D and [Vikhlinin ct al.| 
( [2006 ), who found slightly srnller f gas . W e w ill discuss 
the possible origin of these differences in § 7.2. 



4.4. Entropy profile 
Asi expcctcd for approximately hydrostatic gas , we find 



that Iwe obtain a good fit to the data, with a. model requir- 
ing a monotonically rising entropy (S) profile. We show 
the model profile in Fig [To] (grey shad ed region - ) scale d 
by the "characteristic entropy" K 50 o ( Voit ct aL] gQQ5[ ) , 
and shown as a function of fraction of R500 reached. In 
the inner part of the system, the slope (S ~ R^) is 
slightly steeper than the canonical (3 ~1.1, and the nor- 
malization is significantly enhanced over the "baseline" 
model for gravity -only structure formation simulations 
( [Voit ct "ah 2005), indicating significant entropy injec- 
tion! Above ^40 kpc (^0.07Rsoo), the entropy profile 
flattens significantly, and converges with the (steeper) 
baseline model by R v ir- For maximum flexibility, we al- 
lowed an additional break at large radii, but found that 
it was poorly constrained, and the overall shape of the 
distribution is very flat from ~40 kpc to ~R V j r ; in fact, 
fits omitting the break are able to fit the data just as 
well, with o nly a minimal impact on the recovered gas 
properties (§ 6.9). 

To provide a less model-dependent view of the entropy 
profile, we overlay in Fig [T^ a series of data-points, which 
are directly computed from the deprojected density and 
temperature profiles. These were obtained by emulating 
the behaviour of the "projet" Xspec model, and correct- 
ing for emission projected into the line of sight from out- 



-<r> 



R / R 500 

Fig. 10.— Entropy profile model of RXJ 1159+5531, scaled by its 
characteristic entropy, and shown as a function of R500 ■ The grey 
shaded region is the 1-cr confidence region determined from our 
mass model fit. We overlay deprojected data-points from Chandra 
(triangles) and Suzaku (stars) , which are determined more directly 
from the data, and agree well with the model. We stress the model 
is not fitted to these data. We note that the deprojection proce- 
dure in large spatial bins likely intmduccs non-neglible (unphysical) 
noise into the data-points (see § p.4| ) , and so they should be treated 
with caution. The dotted line indicates the "bnanlyjp" prprlir-t^nna 
from gravitational structure formation (Voit et al. 2005), and we 
find that the "f gaa -corrected" entropy pronle (blue snaded region; 
see text), which agrees well with the baseline predictions out to 
~R v i r . Also shown (yellow region) is the scaled entropy profile of 
the isolated elliptical galaxy NGC 720, illustrating more entropy 
injection in the lower-mass system. 



side the outermost annulus (see Hll). These data agree 
well with the smooth model, giving us confidence in our 
treatment of the data. 



Following Pratt et al. ( 2010 ), we investigated the ef- 
fect of scaling the entropy profile by a simple correction 
factor, S corr = S' o fc s (/ g /0.17) 2/ ' 3 , where S corr is the cor- 
rected entropy profile, and S Q b s is the observed (s caled) 



distri bution. As for more massive galaxy clusters (Pratt 



et al J) , and the isolated elliptical galaxy NGC 720 ( |Hll| ), 
we find that S corr is in much better agreement with the 
baseline entropy model (Fig [To| ). 

4.5. Rosat surface brightness profile 

We next explored whether the model fitted to the 
Chandra and Suzaku data are c onsistent with the Rosat 
surface brightness profile (e.g. Eckert et al. 2011). To 



do this, we computed the three-dimensional gas emis- 
sivity, based on our best-fitting models for the temper- 
ature, density and abundance profiles. This model was 
projected onto the sky and folded through the appro- 
priate Rosat PSPC instrumental responses. To account 
for possible mis-calibration between the satellites, we al- 
lowed an arbitrary scaling of the model normalization 
between Chandra and Rosat. We broadened the surface 
brightness model by folding in the instrumental PSF, 
evaluated at 1 keV. We added a constant (sky) back- 
ground component and allowed its normalization, and 
the aforementioned scaling factor, to fit to the Rosat sur- 
face brightness profile, using dedicated software based 
around the MINUItQ fitting library. The best-fitting 

18 http: / /lcgapp.cern.ch/project / els /work- 

packages / mathlibs / minuit /index. html 
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Fig. 11. — Comparison of the flat-fielded, background-subtracted 
0.3-7.0 keV Chandra (black squares) and 0.42-2 keV Rosat (blue 
triangles) surface brightness profiles for RXJ 1159+5531. Overlaid 
are the profiles predicted by the best-fitting model; the shaded re- 
gions correspond to the l-tr uncertainty. The vertical displacement 
between the Rosat and Chandra data reflects differences in the ef- 
fective area curves and energy-bands. In practice, our model agrees 
very well with the Rosat data, out as far as ~17'. The vertical lines 
indicate R500 and R v i r - 

value of the scaling factor (0.93 ± 0.03) indicates good 
overall agreement, although there may be a modest cal- 
ibration discrepancy, at least when observing a ~1 keV 
source at large (~17') off-axis angles with Rosat. Nev- 
ertheless, such a modest discrepancy will not affect our 
conclusions in the outer part of the group. The surface 
brightness profile was fitted out to ~25', and became 
consistent with the background outside ~12'. 

In Fig [n], we show the background-subtracted Rosat 
surface brightness profile (triangles) and, for comparison, 
the Chandra ACIS-S3 data. We overlay the predicted 
surface brightness models, illustrating excellent agree- 
ment with the Rosat data out beyond R v i r - This strongly 
supports our treatment of the Suzaku (and Chandra) 
data, and confirms that the entropy profile does not ex- 
hibit substantial flattening outside ^Rsoo- 

5. THE COSMIC X-RAY BACKGROUND 

As is clear from Fig g, the accurate determination of 
the gas properties in the outer two Suzaku annuli requires 
the b ackg round to be determined with high accuracy (see 
also § 3.2). The dominant background component of rele- 
vance is actually the cosmic X-ray background (CXB) re- 
sulting from (unrelated) undetected, background point- 
sources. In this section, we explore the accuracy with 
which the CXB component has been fitted in our Suzaku 
analysis. 

Chandra's spatial resolution allows a significant frac- 
tion of the CXB to be resolved into individual point 
sources, at least near the optical axis. By disentangling 
them from any diffuse emission, resolving the sources in 
this way allows the CXB spectral shape and normal- 
ization to be determined more precisely. Deep Chan- 
dra observations have yielded accurate measurements of 
the logN-logS distribution (and hence, average surface 
brightness) of bac kground point sources along particu- 
lar sightlines (e.g. Luo et al. 2008 ). Alternatively, X-ray 
spectra from regions of the sky free of foreground contam- 
ination have allowed the CXB shape and normalization 




Energy (keV) 

Fig. 12. — Composite spectrum of the detected Chandra point- 
sources in the region of interest with fluxes < 2xl0 14 crg s -1 cm -2 . 
The brightest few sources were omitted so that they do not unduly 
bias the fit. Overlaid is the canonical powerlaw (r = 1.41) used to 
fit the CXB component, which gives a reasonably good overall fit 
to the data. 




Fig. 13. — The measured differential log N-log S relation for 
point-sources within the Chandra field of view. The down-turn 
at low fluxes is due to source detection incompleteness; we fold in- 
completeness and Eddington bias corrections into the fitted model. 
The dottqd red line is the expected fit, based on the hard-source 



Luo et al 



(2008), and the solid black line is the CDF-S 



best-fitting model, rescaled to fit the data 

to be careful ly calibrated, averaged ove r small portions of 
the sky (e.g. De Luca & Molendi 2004). However, due to 
cosmic variance and stochastic effects, the inferred sur- 
face brightness will not be well-enough known to use this 
information along other, arbitrary sightlines, such as that 
to RXJ 1159+5531. Ideally, therefore Chandra should be 
used to resolve the point-sources in each of the Suzaku 
annuli, and the resultant spectra can be used to refine 
the Suzaku analysis. Unfortunately, the more obstructed 
field of view of Chandra (at least in the ACIS-S con- 
figuration), coupled with the substantial degradation of 
the PSF (and hence, reduced detection efficiency) at > 4' 
away from the optical axis, means that multiple Chandra 
observations are needed to mosaic the entirety of each 
Suzaku annulus at high enough precision. To compli- 
cate matters further, since point sources can be variable, 
contemporaneous observations with each satellite would 
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ideally be used. Nevertheless, while individual sources 
can vary significantly, the integrated source properties 
are not expected to b e strongly a ffected by variability 
(e.g. |Kraft et~aT| |200l| ; |Zezas et al.||2004| ). 

Since the existing Chandra data did not meet these re- 
quirements, it was not possible to improve our Suzaku 
constraints using this approach. Nevertheless, it is still 
important to verify consistency between the Chandra and 
Suzaku CXB measurement. To do this, we first examined 
the sources detected by Chandra that lay within each 
Suzaku region. In the outermost region, we detected 9 
sources, which we estimate (below) to contribute only 
^20% of the flux within this aperture. This estimate 
reflects both the fact that only ^50% of this region is 
exposed with Chandra, and the large off-axis angles (~9- 
13') under scrutiny. In the 2-5' aperture, however, the 
Chandra data were more helpful; we estimate that ^80% 
of the CXB flux was resolved into 26 detected sources. 
These estimates, of course, assume that the detected 
sources did not vary in brightness significantly between 
the Chandra and Suzaku observations. 

In principle, one can sum the spectra of all detected 
sources within each region, and use that to determine 
the shape (if not the normalization) of the CXB spec- 
trum. However, this estimate suffers from a bias, since 
it assumes that the spectrum of a background AGN is 
independent of its flux, which is not true. If we omit- 
ted sources above a particular flux limit from this calcu- 
lation, the composite spectrum was systematically flat- 
ter. Fits to the CXB spectrum measured from unre- 
solved background emission in source-free regions of the 
sky should not be affected by this problem, and so in 
our default analysis, we parameterized the CXB spec- 
trum as a powe r law w ith T = 1.41, as found by De 
Luca fc Molendj ( 2004 ). The incompleteness correction 
of the composite spectrum that is necessary to test formal 
consistency between this model and the Chandra data 
is beyond the scope of this paper. However, in Fig |l^ 
we show the composite spectrum of sources fainter than 
2 x 10~ 14 erg s _1 cm~ 2 (omitting the few, bright point 
sources in the field so as not to skew unfairly the spectral 
shape), which agrees r easo nably well (if not perfectly) 



with this model. In S 6.2 we find that modest varia- 



tion in the slope of the CXB component (±5%) does not 
strongly affect our conclusions. 

While the CXB shape may vary modestly between 
sightlines, the uncertainty on its normalization is likely 
to be more problematical for our analysis. To verify con- 
sistency between the Suzaku flux in each aperture and 
the Chandra esimates requires an estimate not only of 
the mean flux, but also the uncertainty on it, from un- 
detected sources in each region. To measure the mean 
flux, it was first necessary to determine the normaliza- 
tion and shape of the average logN-logS distribution for 
the point sources along the RXJ 1159-1-5531 line of sight. 
Since cosmic variance is unlikely to be strong between 
each Suzaku region, we considered all sources detected by 
Chandra t hat lay within the S u zaku X ISO field of view. 
Following Humphrey & Buote (2008), we subtracted a 
local background from each source spectrum and con- 
verted the counts to 0.5-7.0 keV flux, assuming a power- 
law model (r=1.55) with Galactic absorption, while tak- 
ing into account the spatial dependence of the effective 
area. Given the measured fluxes and spatial distribu- 
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Fig. 14. — The 0.5-7.0 keV surface brightness profile (Ecxb) of 
the CXB component, measured from the Suzaku, Rosat and Chan- 
dra data. We show the best-fitting T,cx B from the default Suzaku 
analysis (black dash-dot-dot line), and Ticxb determined indepen- 
dently from the Suzaku data in each region (black crosses). Also 
shown are the estimates from the Rosat surface brightness profile 
(dashed error bars; blue circles) and from fitting the Chandra point 
source distributions (red triangles). The blue dotted line is the ex- 
pected flux from integrating theJogiiLkigS point source distribution 



along the CDF-S line-of-sight (Luo et al 
excellent agreement between the different 



2004 see text). We find 
estimates. Note that the 



CXB component need not be constant with radius, and so we only 
expect good agreement only between the fitted data-points, not 
the approximate (constant) models shown. 

tion, we do not believe any of these sources to be X-ray 
binaries associated with the central galaxy. The result- 
ing logN-logS distribution (in dN/dS form) is shown in 
Fig [UJ In order to interpret these data, it was necessary 
to fit a model, correcting for both i ncompleteness at low 
fluxes and the Eddington bias. In Humphrey & Buotc 



d2008|) 



we discuss in detail how these corrections were 
carried out. In short, we corrected the model {not the 
data) based on the results of extensive Monte Carlo sim- 
ulations in which fake sources were added (in a spatially 
uniform, random fashion) to the Chandra images, and 
the source detection algorithm was used to try to detect 
them and determine their flux. 

We found that the dN/dS distribution could be well- 
fitted with a (incompleteness-corrected) broken powerlaw 
model (with a break at 7.4 x 10~ 15 erg s _1 cm~ 2 , and 
negative logarithmic differential slopes 1.3 and 2.6 below 
and above the break, respectively) that we also found 
co uld fit the h ard-b and logN-logS distribution reported 
by [Luo et alj ( |2008| ) for the Chandra Deep Field South 
(CDF-S) line of sightj^]. The best-fitting model is shown 
in Fig H3, which is very close to the fit to the CDF-S data 
(shown, correctly normalized for the observed region of 
sky). This good agreement gives us confidence in our 
treatment of the data. 

While the average normalization of the CXB compo- 
nent along the RXJ 1159+5531 sightline can be deter- 
mined in this way, the actual flux measured in any annu- 
lus is subject to stochastic scatter about this value due 
to the discrete source nature of the CXB. To account for 
this, we undertook Monte Carlo simulations, as outlined 
below: 

19 Taking into account the different energy-bands used, and the 
spectral models used to convert counts to flux. 
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1. Using the best- fitting logN-logS relation, described 
above, we determined the total number of sources per 
square degree, No with flux > 10 _16 erg s _1 cm~ 2 . On 
each simulation, we added Gaussian noise to No to reflect 
uncertainties (±12%) in the fit to the logN-logS relation. 

2. For each Suzaku region i, which has area A4, we drew 
Nj sources that have fluxes distributed like the logN-logS 
relation between 10 -16 and 10 -12 erg s _1 cm~ 2 . Nj was 
itself Poisson distributed around NqAi. 

3. Based on the fraction of the Suzaku region covered by 
Chandra, and our source detection incompleteness esti- 
mates at each flux, every source has a particular proba- 
bility of being detected. Using this information, we ran- 
domly assigned a state (detected, or not detected) to 
each source. 

4. In each region i, we estimated the CXB flux by 
summing up the flux of all simulated sources flagged 
as "undetected" , and added in the flux of the real de- 
tected sources and the expected flux of sources below 
10~ 16 crg s _1 cm~ 2 (by integrating the logN-logS rela- 
tion). 

5. We repeated steps 1-4 a large number of times to de- 
termine the average CXB flux and its uncertainty in each 
region i. This estimate accounts for both stochastic scat- 
ter and the uncertainty on the logN-logS normalization 
from our fit. 

In Fig |l4], we show the estimated X-ray flux for each bin. 
along with the value measured from the Suzaku spec- 
tral fitting. We show both the best-fitting value for our 
default analysis (in which the CXB normalization per 
square degree is tied between all the annuli), and the 
constraints on the CXB flux when the the normalization 
is allowed to fit freely in each region. The Chandra re- 
sults appear to be in good agreement with the Suzaku 
measurements^. If the integrated flux is very sensitive 
to the properties of individual, bright point sources, we 
would expect the CXB surface brightness measured by 
Suzaku to show significant scatter between each radial 
bin, in excess of the measured statistical error. In fact, 
the CXB surface brightness is consistent with being con- 
stant with radius. 

As a final consistency check, we also explored whether 
the Rosat data agreed with the Suzaku results. We ex- 
tracted the Rosat radial surface brightness profile in the 
band 1-2 keV; based on the best-fitting Suzaku back- 
ground model, we expect that ^95% of the background 
emission in this region comes from the unresolved CXB. 
We restricted the photons to be from the region of sky 
imaged by XISO (although we also excluded a 3.4'-wide 
region in the vicinity of the inner support ring shadow). 
We subtracted off the particle background, the expected 
(small) contribution from the soft Galactic background 
components (which was assumed to be spatially uni- 
form), and the surface brightness pr ofile predicted from 
our best fitting mass model (see § 4.5). The remain- 
ing counts should come from the CXB emission. Rebin- 
ning the profile to match the Suzaku annuli, and convert- 
ing counts to flux by folding the canonical CXB model 
through the Rosat response matrices (computed close to 

20 We note that the slight difference in the powerlaw slopes used 
to flux the sources (T = 1.55) and to fit the CXB (r = 1.41) leads 
the Chandra data to underestimate the Suzaku flux only by ~6%, 
which is far smaller than the statistical errors. 



the centre of RXJ 1159+5531), we obtained an estimate 
of the CXB surface brightness in excellent agreement 
with the Suzaku measurements (Fig |l4| ) . 

6. SYSTEMATIC ERROR BUDGET 

In this section, we address the sensitivity of our re- 
sults to various data analysis choices that were made. In 
most cases, it is difficult or impractical to express these 
assumptions through a single additional model param- 
eter over which one might hope to marginalize, and so 
we adopted the pragmatic approach of investigating how 
our results changed if the assumptions were arbitrarily 
adjusted. We focused on those systematic effects likely 
to have the greatest impact on our conclusions, and list 
in Tables [l] and || the change to the marginalized value 
of each key parameter for each test. We outline below 
how each test was performed. 

6.1. Dark Matter profile 

One of the major sources of uncertainty on the re- 
covered mass model is the coice of DM mass model. 
While the NFW model is theoretically motivated, we 
also experiemente d with the so-called "co red logarith- 
mic" mass model ( Binney fc Tremaine 2008 ). This model 
tends to predict higher masses (by ~60%), and corre- 
spondingly smaller gas fractions, at large scales. Al- 
though we cannot distinguish between the NFW and 
cored logarithmic models on the basis of % 2 alone, the 
ratio of the Bayesian evidence (2.6 x 10 -3 ) implies that 
the cored logarithmic model, with the adopted priors (a 
flat prior on the asymptotic circular velocity, between 10 
and 2000 km s _1 , and a flat prior on logior c , where r c 
is the core radius, over the range < logior c < 3.) is a 
poorer description of the data at ~3.0-o\ 

Another modification to the DM halo profile that is 
well-motiva ted theoretically, although less secure qbser 



vationallv f|Humphrev et al. 2006; Gnedin et al. 2007; 



Napolitano et al! |2010|). is so - called "adiabatic contrac- 



tion" (jBlumcnthai ct aL||l986| ; jGnedin et al.||2004| ; [Abadi 
;t al. 2009|) , where the DM halo density profile reacts to 



the gravitational influence of baryons that are condens- 
ing into stars by becoming c uspicr. Modify i ng th e NFW 
profile with the algorithm of |Gncdin ct al. ( 2004j )p], has 
only a very slight effect on the best-fitting mass model 
("AAC" in Tables and |). This reflects that the scale 
radius of the DM halo is much larger than the effective 
radius of the stellar mass component. 

6.2. Background 

Since the data were background-dominated in the out- 
ermost Suzaku annuli, the treatment of the background 
was a potentially serious source of systematic uncer- 
tainty. To investigate the extent to which our results are 
sensitive to this, we explored a range of different choices. 
First, for the Chandra data, we adopted the standard 
blank-field events files distributed with the CALDB to 
extract a background spectrum for each annulus. Since 
the blank-field files for each CCD have different expo- 
sures, spectra were accumulated for each CCD individu- 
ally, scaled to a common exposure time and then added. 

21 Using the CONTRA code publicly available from 
http: / /www. astro, lsa.umich.edu/~ognedin/contra/ 



16 



The spectra were renormalized to match the observed 
count-rate in the 9-12 keV band. These "template" spec- 
tra were then used as a background in Xspec, and the 
background model components were omitted from our fit. 
This gave a formally poorer fit, but did not strongly af- 
fect our conclusions. Second, since the Suzaku non X-ray 
background spectra generated by xisnxbgen may be un- 
certain, we allowed their normalization to scale by ±5%. 

At the temperature of the gas in the outermost annuli 
(~1 keV), the dominant background component is the 
CXB. In §|, we demonstrated that the best-fitting CXB 
model to our data is in excellent agreement with predic- 
tions for the line of sight to RXJ 1159+5531. However, 
there were still uncertainties in our treatment. Specifi- 
cally, by default, we assumed a constant surface bright- 
ness for the CXB component, which may not be formally 
correct. We therefore experimented with allowing the 
CXB normalization to fit freely in each Suzaku annulus 
(Fig [li]). This did not significantly affect our results. 
Additionally, there is uncertainty on the spectral shape 
of the CXB component. Given the statistical uncertainty 
on t he shape of t he CX B parameterization found by Dc 
Luca | fc Molendj ( ^004j ), we varied the slope of the CXB 
powerlaw component by ±5%. Although this affected 
the error on the gas density at the largest radii, our con- 
clusions were not significantly altered. Allowing a larger 
change in the slope did have a more significant effect on 
the density and temperature, however, motivating the 
need for deep Chandra data to res olve and constrain the 
CXB at the largest scales (§ |7.6| ). We summarize the 
results in Tables]!] and || ( "ABackground" ) . 

6.2.1. Solar Wind Charge Exchange 

An additional background component can arise from 
the interaction of the Solar wind with interstellar ma- 
terial and the Earth's exosphere. This should manifest 
itself as a time- variable, soft component that can be mod- 
elled as a series of narrow Gaussian lines, the inten sity of 



which correlate with the Solar wind activity (e.g. 3now- 
' en ct al | |2004|) 



We did not explicitly include compo- 
nents to account for this "Solar Wind Charge Exchange" 
(SWCX) in our fits, although the soft background com- 
ponents are partially degenerate with it. To explore 
whether the SWCX could have affected our conclusions, 
we used the Solar Wind Ion Composition Spectrometer 
(SWICS) instrument aboard the Advanced Compo sition 



to 



Explorer (ACE) spacecraft ( McComas et al. 1998 ) 
identify periods of enh anced SWCX emission. Following 
[Bnowdcn ct al] ( [2004 ), we assumed it to be negligible 
when the +7 /0 +B ratio falls below ^0.2, and selected 
times during each observation which met this criterion. 
While the Chandra data were moderately affected (~60% 
of the data were excluded by this cut), in the low sur- 
face brightness regime, the Suzaku data were most im- 
portant, and these were only mildly affected (~9% of the 
data were removed by this cut; the 0.5-1.0 keV count- 
rate was enhanced only by ^6% if all the data were used). 
The corrected Chandra and Suzaku spectra were fitted to 
obtain the temperature and density profiles, and folded 
through our mass modelling apparatus. We found that 
our results were not substantially affected by correcting 

22 Based on the publicly released data available from 
http://www.srl.caltech.edu/ACE/ASC/index.html 



for the SWCX component ("ASWCX" in Tables | and 

6.3. Stray light 

Stray light from bright point sources within ~ 1° can, 
in principle, be scattered by the mirrors into the field of 
view, and provide an additional source of background. 
By far the brightest point source identified by the Rosat 
All-Sky Survey (RASS) within this angular distance from 
RXJ 1159+5531 is is NGC 3998, a foreground Syl galaxy 
which is 17' from RXJ 1159+5531. To investigate the po- 
tential impact of stray light, we included an additional 
component in our spectral modelling, corresponding to 
the contamination from NGC 3998 to each annulus. We 
estimated the amount of stray light leakage by using 
xissimarf gen to generate suitable ancillary response 
files for each annulus, with the locus of NGC 3998 as 
the origin for the photons. We modelled the emission 
from NGC 3998 as an absorbed powerlaw, with Njj, T 
and the normalization fixed to the values derived from 
archival XMM-Newton data by |Ptak ct al] fl2004j). We 



found that the stray light contamination was < 20% of 
the source flux below < 2 keV. The stray light compo- 
nent is much harder than the ~1 keV spectrum of the 
hot gas, and is partially degenerate with the CXB spec- 
tral component. Therefore, we found that adding this 
component had minimal effect on our results ("AStray 
light" in Tables0 and |). 

6.4. Projection/ Deprojection 

In our analysis, we modelled the projected tempera- 
ture and density in each annulus by evaluating the hy- 
drostatic model for the temperature and density in three 
dimensions, and projecting it onto the line of sight. In 
principle, the results can be sensitive to the outermost 
radius used in the projection calculation. By default, we 
adopted 2R v i r , but we also explored varying this limit 
between R v j r and 3R v h . This had a minor impact on our 
results ("Ai? ma:r " in Tables |] and |) . 

In this work, we fitted the projected, r ather than the 
deprojected data (as done, for example, in H06 ) . In gen- 
eral, fitting the projected data leads to smaller statistical 
error bars, but potentially larg er systematic uncertainties 
(e.g. Gastaldcllo ct al. 2007b), and so it is important to 
investigate the likely magnitude of such errors. To do 
this, we examined the effect on our results of spherically 
deprojecting the data. We achieved this by using the 
Xspec projet model^J. To account for emission projected 
into the line of sight from regions outside the outermost 
annulus, we added an apec plasma model to each annu- 
lus, with abundance 0.2 (consistent with the outermost 
annulus) and the temperature and normalization deter- 
mined from projecting onto the line of sight the best- 
fitting gas temperature and density models described in 



4.1 



, but considering the models only outside ^1200 kpc. 
We note that the results from a deprojection analy- 
sis in the large radial bins used here should be treated 
with considerable caution. Since the procedure assumes 

23 In practice, it was more convenient to emulate the behaviour 
of projet by adding multiple "vapec" plasma models in each annn-. 
lug with the relative normalizations tied appropriately (e.g. Kriss 
pt al] |l983| ) . This allowed data from the multiple Suzaku instru- 
ments to be fitted simultaneously. 
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constant density and emissivity (hence temperature and 
abundance) in each shell, which is substantial simplica- 
tion (see Fig|J), this ca n introduce significant, unphysical 
nois e (e.g. see § 3.3 of Buotc 2000a ; Finoguenov fc Pon- 
i |1999| ). Nevertheless, wc found that the best- fitting 



derived results were not significantly affected when the 
deprojected (rather than the projected) temperature and 
density profiles were used ( "ADeprojection" in Tables |l| 
and P). In Fig |, we show a series of mass "data-points" 
obtained from the deprojected density and temperature 
profiles, usi ng the "traditional" ma ss analy sis method 



desc ribed in Humphrey et al. ( 2009a , see also Humphrey 
& Bj^ti] |20iq ). These data agree very well with the 
best-fitting mass model found in our projected analysis. 
Similarly, the entropy profile derived directly from the 
deprojected data displays overall agreement with the re- 
sults of our projected analysis, although the data-points 
exhibit unphysical "choppiness" due to the deprojection 
noise (Fig|l(i|). 

6.5. Priors 

Since the choice of priors on the various parameters is 
arbitrary in our analysis, it is important to determine to 
what extent they could affect our conclusions. To do this, 
we replaced each arbitrary choice in turn with an alter- 
native, reasonable prior. Specifically, for each parameter 
describing the entropy profile, we switched from a flat 
prior on that parameter to a flat prior on its logarithm. 
We used a flat prior on the DM halo mass, rather than on 
its logarithm, and, instead of the flat prior on log cdm, 
we a dopted the di s tribut ion of c around M f ound by ei- 



ther |Buote et a"j] ( [2007D or |Macci6 et all ( |2008[ ) as a 
(Gaussian) prior. The effect of these choices is no larger 
than the statistical errors on each parameter, especially 
for the baryon fraction measured at i?2oo or higher over- 
densities ("AFit priors" in Tables @ and §). 

6.6. Stellar light 

Since the scale radius of the DM halo is very much 
larger than the effective radius of the stellar light, we 
would not expect a careful treatment of the stellar light 
to be important for accurately measuring the total mass 
of the system ( in contrast with the lo wer-mass, elliptical 
galaxy, regime: Humphrey et al. 2006 ). To illustrate this, 
we excluded the central ^20 kpc (roughly twice R e ) from 
our analysis. Since the stellar mass is relatively unimpor- 
tant at this radial scale (Fig@), we fixed the stellar M/L 
ratio to its best-fitting value. This had a minimal impact 
on our conclusions. 

Our measurement of fb includes a canonical amount 
of intra-cluster light, which was not directly detected. 
Since the true contribution of this component (and its 
M/L ratio) is unknown, we explored how significantly fb 
was affected if this component was completely omitted, 
or if its M/L ratio was fixed at 1 (slightly higher than 
the best- fitting M/L ratio for the central galaxy). These 
choices did not substantially affect our results ( "AStars" 
in Tables and |). 

6.7. Emissivity correction 

In our default analysis, the projected temperature and 
density profile were weighted by the gas emissivity, folded 
through the instrumental responses (for details, see Ap- 



pendix B of ^astaHcIIocFar 2007b). Since the computa- 
tion of the gas emissivity assumes that the three dimen- 
sional gas abundance profile is identical to the projected 
profile (which is unlikely to be true), it is important to 
assess how sensitive our conclusions are to the emissivity 
correction. To do this, we adopted the extreme approach 
of ignoring the spatial variation of the gas emissivity al- 
together. We found that this had a very small effect on 
our results (AWeighting in Tables [l] and ||). 

6.8. PSF correction 

Our Suzaku analysis results depend on the careful 
treatment of spectral mixing between each annulus, 
which occurs due to the modest spat i al re solution of 
the telescope (see Ell; Rciprich et al. 20091). Our ap- 
proach was to calculate mixing between the annuli with 
the xissimarf gen task, which uses ray-tracing. To ex- 
plore whether small errors in this procedure could affect 
our results, we experimented with scaling the amount of 
light that is scattered into each annulus by ±5%. This 
did not appreciably affect our conclusions ("APSF" in 
Tables |l]and §. 

6.9. Remaining tests 

We here outline the remaining tests we carried out, 
as summarized in Tables [j] and ^. First of all, since 
the inter-calibration of the Suzaku XIS units may not 
be perfect, we experimented with using only one of the 
units in the Suzaku analysis, and cycled through each 
choice. We also considered using only the Chandra data. 
We found that our results were resilient to this choice 
("AInstrument"). 

To assess the impact of various spectral-fitting data 
analysis choices on our results, in turn we varied the neu- 
tral hydrogen column density by ±25%, performed the 
fit over a restricted energy range (0.7-7.0 keV), and re- 
placed the APEC plasma model with a MEKAL model. 
The impact of these choices is comparable to the statisti- 
cal errors on the parameters ("ASpectral"). To examine 
the error associated with distance uncertainties, we var- 
ied the distance to RXJ 1159+5531 by ±30%, finding the 
effect, particularly on fb and f gas to be relatively minor 
("ADistance"). 

In order to provide maximum flexibility in our fitting at 
large radius, our parameterization of the entropy profile 
included an arbitrary break at large radius. We found 
that, in general, the break radius tended to large scales, 
consistent with it not being required in our fits. We 
experimented with fitting a model without this break, 
obtaining a good fit (x 2 /dof=14.6/17) in close agreement 
with our default model ("AEntropy"). 

Finally, to examine the possible errors associated with 
our treatment of the covariance between the density 
data-points, we investigated adopting a more complete 
treatment that considers the covariance between all the 
temperature and density data-points, as well as adopting 
the more standard (but incorrect) approach of ignoring 
the covariance altogether. The best-fitting \ 2 value did 
not change significantly (±5) with these different choices, 
and so the effect on our fits is not large ( "ACovariance" 
in Tables and |). 

7. DISCUSSION 
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By combining our new deep Suzaku observation of the 
very relaxed fossil group/ poor cluster RXJ 1159+5531 
with archival Chandra data, we have now obtained an 
unprecedented census of the baryons and dark matter 
over almost 3 orders of magnitude in radial range, from 
~kpc scales to the virial radius (-Rios = 1100 kpc). We 
discuss here the implications of these new measurements, 
in detail. 

7.1. Hydrostatic equilibrium 

Our best-fitting hydrostatic model fits the density and 
temperature data-points extremely well, strongly sup- 
porting the hydrostatic approximation for this system; 
despite highly nontrivial temperature and density pro- 
files, a smooth, physical mass model and a monotonically 
rising entropy profile (required for stability against con- 
vection) were able to reproduce them well. This would 
require a remarkable conspiracy between the tempera- 
ture and density if this approximation were seriously in 
error. At the crucial, largest scales, it is striking that 
we see no evidence of the peculiarities, such as a flat en- 
tropy profile or very high f gas value, that c haracterized 
the distur b ed systems Perse us and Virgo (Bimionescu 
\it al. 2011; Urban et al. 2011), and might point to devia- 
tions from the hydrostatic approximation in those sys- 
tems. This agreement between the models and data 
holds even down to ~kpc sca les, which had been omit- 
ted iii some previous studies (IVikhlinin ct all 1 20061; Bun| 



dist urbance within the central ^10 kpc (tj 2.1). recent 



studies of systems that are manifestly more 



isturbed 



than RXJ 1159+5531 indicate that the mass derived from 
azimuthally-averaged temperature and density profiles 
is relatively r obust to such features (global errors be- 
ing ~1 0-20%: |Churazov et al.||2008| ; [Buotc fe Humphrey 
|2011a, for a detailed discussion). In any case, whether 



or not we include data at these sm all s cales has only a 
minimal effect on our conclusions (§ 6.6). 

The closeness of the system to hydrostatic equilibrium 
is unsurprising, given its round, relaxed-looking X-ray 
isophotes (Fig [j]). Numerical simulations of structure 
form ation suggest that deviations from hydrostatic eaui- 



librh m in morphologically relaxed-looking systems are 



tion should be no larger than ^25% (e 



Buotc 
larnim 



Tsa 
2TO 



i | HH |Nagai et a l 



Tsai et al. 


1994; 


Fittarctti & Val- 



''inally, it is important 



ang et al, 

to note that, as a (highly evolved) fossil group, the accre- 
tion rate of the system should be small, hence possible 
associated morphological and dynamical perturbations 
at the largest scales will be minimized. 

7.2. Mass Profile 

The total mass profile of RXJ 1159+5531 out to the 
virial radius is very well fitted by a model comprising 
an NFW dark matter halo, a stellar component and gas 
(plus, at the smallest scales, a black hole), in good ac- 
cord with theory. We note that such agreement is not 
guaranteed for any model, since at least one other model 
(the cored lo gari thmic potential) is rejected at ~3-<r in 
our fitting (§ 6.1). Previous studies have suggested that 
fossil groups may have significantly enhanced concentra- 
tion paramet ers, likely due to their ea rly epoch of for- 
mation (e.g. Khosroshahi et al. 2007 ). In the case of 



RXJ 1159+5531, however, the halo concentration param- 
eter (11.2 ± 1.6) is in good agreement both with theoret- 
ical predictions ( from dark matter-only simulations) for 
this mass range ( Bullock et al. 2001 ; Maccio ct al. _ 2008|) 



and w ith the empirical relation derived by 



( |2007| F. Similar results were obtained for t 



Buote et al 
re other fos- 



sil systems discussed in G07 and 3un ct al. (2009). 

Our best-fitting mass model parameters are m excel- 
lent agreement with the Chandra measurements made by 
G07 , although with much smaller error bars (Fig@). The 
reduction in the size of the errors reflects both the addi- 
tion of the high-quality Suzaku data at large radii (which 
helps to pin down the flattening of the mass profile at 
these scales), and various recent enhancements made to 
our mass-modelling procedure that improve the robust- 
ness of the fit to sources o f systematic uncertainty (see 
Buote & Humphrey 2011a). In particular, by rigorously 
enforcing the Schwarzschild criterion for stability against 
convection, the parameter space is su bstantially shrunk. 

While our results agree with G07, puzzling discrep- 
ancies e xist with the other published mass par ameteri- 
zations QVikhlinin ct al.| |2006| ; pun ct al.|[2009| ). These 
fits used only the Chandra data, and eliminated the 
central ^40 kpc from the fit. The mass distribution 
was inferred from the gas density and temperature pro- 
files using the traditional " smoothed inversion" method 



|ot al 2009|). While there is weak evidence asmali et all |2009a| 



( Buotc fe Humphrey! pO~lla| , for a review), w hich can in- 
troduce non-negligible systcmat 



In Fig 



/stematic errors (e.g. Humphrey 
we show the (projec ted) tem- 
perature and (deprojected) density profiles of Vikhlinin 
ct al.| ( |2006j) , while we show their recovered mass pro- 
file in Fig |[ While the temperature and density data 
agree well with our results, there is modest disagree- 
ment in the mass distribution at large (> 300 kpc) and 
small (< 50 kpc) scales. We attribute these discrepan- 
cies to systematic errors in the smoothed inversion tech- 
nique, but they cannot explain the significant differences 
in the parameterized fits. Fitting an NFW model to these 
data, we obtain best-fitting parameters (log A/500 = 13.8, 
logesoo = 0.8) that are very close to our measurements, 
albeit with large fit residuals at large scales. Ra ther than 
fitting the enclosed mass distribution directly, Vikhlinin 



et al. actually fitted the gravitating mass density profile 
derived from their data. In this case, an NFW density 
model fit was primilary constrained by the data at large 
radii, yielding a much lower concentration. In both cases, 
however, a simple powerlaw actually fitted the profiles 
better than NFW. We therefore attri bute the lower c on- 
centration and higher mass found by Vikhlinin et al. to 
systematic errors in their derivation of the density pro- 
file, which is inherently more uncertain than the enclosed 



mass. Since sun et al. did not provide individual mass 
profiles in their paper, it was not possible to carry out a 
similar comparison with their work. 

7.3. The baryon fraction 

We were able to place tight constraints on both the 
gas and baryon fraction out to the virial radius of the 

24 Strictly speaking, this comparison is not independent, since 
R V 1 1 1 ^a-|-5^i w as one of the systems used to determine the 



Buote et al. (2007) relation, albeit with much poorer data. Nev- 
ertheless, itXj 1159+5531 does not appear to be unusual in com- 
parison to the other systems used by these authors, and so it is 
unlikely that this relation is largely driven by this one data-point. 
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R/R Iir 

Fig. 15.— The radial profile of f b and f gas for RXJ 1159+5531 
(dark grey, dark blue regions, respectively), versus fraction of R v i r - 
In comparison, we also show the same quantities for the isolated el- 
liptical galaxy NGC720 (light grey, light blue, respectively). Note 
the overall similarity in the f^ profile shape, despite significant dis- 
crepancies in fgas- 

system. While f gas rises fairly steeply with radius (so 
that ^60% of the total gas mass lies outside ^R.500), we 
found that fb rises more modestly from ~0. 1 at ^100 kpc 



to ap p roach the Universal b aryon fraction (Dunkley et al 



|2009| ; |Komatsu et al.| |201l| ) by R vir . This is consistent 
with the picture that X-ray bright gro ups may be bary- 



onically closed ( Mathews et al.| 2005). Intriguingly, in 
the isolated, Milky Way-mass elliptical galaxy NGC 720 
(almost ~2 orders of magnitude lower in mass), we found 
a remarkably similar trend. In Fig [15], we compare the 
profiles for the two systems, scaled with respect to the 
virial radius. Despite the different f ; a S profiles, fi, i s quite 
consistent over a wide radial range. Giodini ct al. ( 2009] ) 
similarly found a conspiracy between the stellar and gas 
mass to produce a total baryon fraction at R 50 o that 
is relatively insensitive to the virial mass of the halo. 
This indicates that, whatever processes are involved in 
redistributing the baryons in the system, in general the 
baryons are not being ejected, at least down to ~Milky 
Way masses. 

While fb asymptotes to the Universal value by R v i r , 
the local gas fraction exceeds this limit outside ^Rsoo- 
Such an enhancement is to be expected if feedback is re- 
distributing baryons in the potential well, but not com- 
pletely evacuating them from the system. The apparent 
delicate balance between the gravitational potential and 
the energy injection from feedback that allows, over a 
wide range of virial masses, the gas to be pushed out 
almost to Rvir, but not evacuated from the system, may 
simply reflect that the mass of gas available to fuel feed- 
back scales with the mass of the system. 

Within R2500 and R.^nn. our f ffaa measurements were in 




100 

Radius (kpc) 

Fig. 16.— Cooling time profile of RXJ 1159+5531. Most of the 
gas in the system has a cooling time longer than the Hubble time 
(indicated with the dashed line). 

we found no evidence of a systematically over estimated 
f gas , which might hint at additional processes such as 
deviations from hydrostat ic equilbrium or spheri city, or 
dumpiness of the IGrM ( gimionescu et al. 201 1| ). The 
low accretion rate expected for an old fossil group may 
reduce the impact of these effects in RXJ 1159+5531. We 
note that we cannot entirely rule out clumping in the 
outer parts of RXJ 1159+5531, but if present, it must 
operate in a finely-balanced conspiracy with fb in such a 
manner as to produce a total baryon fraction very close 
to the cosmological value. In our opinion, such a scenario 
falls foul of Occam's razor. 

7.4. Entropy profile 

The entropy profile of RXJ 1159+5531 shows remark- 
ably s imilar characteristi c s to th a t of other galaxy groups 



(e.g. Gastaldcllo ct al.l 2007a 




Mahdavi et al. 


2005; 


Finoguenov ct al.| 2007; pun ct al. 


|2009|; (Javagnolo et al. 


2009|; |Johnson ct al 2009; |Flohic et al.| 2011). In the 



r 11 , but offset u pwards from the baseline entropy model, 
Vbit ct al. 2005), then flattens outside ~0.06Rsoo (where 
the cooling time ~ 10 10 Gyr: Fig O), and converges with 
the baseline model by ~Rsoo- As expected for hydro- 
static gas, S rises monotonically with radius, and we 
found little evidence of any further breaks in the en- 
tropy profile. This overall shape is similar to what is pre- 
dicted for models in wh ich both stellar winds and AGN 
contribute to feedback ([McCarthy et al.||2010|) . Signifi- 
cantly, we did not find the characteristic flat tening out- 
side ^Rnnn reported for Perse us and Virgo ( Bimionescu 



good agreement with|G07, although slig htly higher than 
t hose found by jyikhlinin et aL ( [2006 ) and Sun ct al 
( 2009 ) (see Fig ||). This reflects a slightly lower gas den 
sity found by those authors (e.g. Fig ||), at least within 
R-2500- These results are also in good accord with simi- 
lar measurements made in other systems (see Fig 11 of isolated elliptical galaxy NGC 720 (|Hll|), when we scaled 



t al. 2011; Urban et al. 2011), either in the deprojected 
data, or m our fits with a parameterized model (which 
had sufficient parameters to capture such behaviour). 
Therefore, whatever processes are involved in producing 
such a feature cannot be ubiquitous in massive groups or 

clusters (or, at least, fossil groups). 

As seen in massive galaxy clusters flPratt et al. 2010] ) . 



the central parts of groups ( Flohic ct al. pOlT ) and the 



Hll| ), indicating that RXJ 1159+5531, despite being a 
fossil group, does not experience a highly unusual evolu- 
tion of its IGrM. In contrast to the outskirts of Perseus, 



the entropy profile by (/ gas /0.17) 2 / 3 , we found it was 
brought into fairly close agreement with the baseline, 
adiabatic model (Fig [h]). For RXJ 1159+5531, we have 
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Radius/R 5[ 

Fig. 17.— Pressure profile foi RVT11SQ+.^1 



line is the canonical model of Arnaud et al 



The dashed blue 



(201C) 



now confirmed that this relation holds even out to the 
virial radius. Taken together with the local f gas con- 
straints, it seems that the primary effect of feedback on 
the baryons is to redistribute the hot gas within the halo. 
rather t han t o raise its temperature, or eject it (Mathews 

fc Gi^poTTl). 



7.5. The pressure profile of RXJ 1159+5531 

Recently, it has become evident that the thermal pres- 
sure profiles of galaxy clusters and groups, when scaled 
appropriately, exhibit remarkable uniformity, at least be- 
tween ~0.5-1.0R,5oo, which has important implications 
for un derstanding and cal ibrating standard scaling rela- 
tions. Arnaud et al. ( 2010) proposed a canonical pressure 
profile model, based on observed clusters wit hin ~RRon 



and n umerical simulations at larger scales. Sun et al 
( 2011 ) confirmed that this profile holds for galaxy groups 
within ^Rsoo. With RXJ 1159+5531, it is now possible 
to test this model at larger scales. In Fig 17, we show 
the scaled pressure profile of RXJ 1159+5531. with the 



canonical model of |Arnaud et al 
scales (<0.1R5oo) 



overlaid. At small 
there is some disagreement, but this 



is not unexpected, based on the cluster and group sam- 
ples. There is general agreement between ~0.2-0.5R5oo, 
but the model significantly diverges from the observed 
pressure distribution at larger scales. Since the global 
properties^ of RXJ 1159+5531 are in agreem ent with 



standard scaling relations (e.g. Sun et al. 2009), it is un- 
likely that this represents an artefact of, for example, its 
fossil group nature. More measurements of the pressure 
profile out to ^R v i r are clearly needed. 

7.6. Constraining the CXB 

Despite the data being background dominated at large 
radius, our results appear to be reasonably robust. This 
is implied by the well-behaved nature of the density, 
temperature and mass profiles, the good agreement be- 
tween the different satellites (in particular, the agree- 
ment with the Rosat surface brightness profile) , and the 
reliable measurement of the CXB normalization in our 

25 M 50 o = (5.9±0.5) X 1O 12 M ; fcT 50 o = 1.78 + 0.05, where kT 500 
is the emission-weighted temperature within the projected aperture 
(O.l-l)Rsqo; Y x = (9.1 ±0.6) x m" Af c h P v whore Y x is the gas 
mass within R500 times kT^oo; see 3un et al. (20091). 



Suzaku fits. Nevertheless, the gas density at the largest 
scales was sensitive to the precise background, specifi- 
cally the CXB spectrum (which is the dominant back- 
ground component). There are, therefore, real concerns 
that, if the adopted model is not accurate, we may be 
mis-characterizing the gas properties in the outskirts of 
the group. As discussed in § |^, one way to mitigate 
these worries is to employ a Chandra mosaic of the outer 
Suzaku annuli, enabling us to resolve a large fraction of 
the CXB component into point sources, from which the 
CXB shape and normalization can be directly measured. 
Combining Chandra and Suzaku in this way exploits their 
complementary capabilities (the high spatial resolution 
of Chandra; the low instrumental background and good 
spectral resolution of Suzaku), and should allow the ro- 
bust detection of gas out to R v i r to become routine. 

8. CONCLUSIONS 

In our joint Suzaku, Chandra and Rosat analysis of 
RXJ 1159+5531, we found: 

1. By choosing objects with an advantageous com- 
bination of distance, mass and surface brightness, it is 
possible to detect the gas to R v i r in a single Suzaku 
pointing. This approach provides better observational 
efficiency than the multiple pointings needed to mosaic 
nearby systems, while providing better azimuthal cover- 
age. 

2. Applying this approach to the morphologically re- 
laxed fossil group RXJ 1159+5531, we were able to mea- 
sure the distribution of gas and gravitating matter out 
to Rvir for the first time in a system with a mass as low 
as 1O 14 M . 

3. Within R V i r , the total baryon fraction approached 
the cosmological value, implying that feedback only re- 
distributed the baryons within the halo, rather than 
ejecting them. 

4. The entropy profile was enhanced over the model 
from gravity-only cosmological simulations, implying sig- 
nificant feedback. However, outside ~0.06R5oo, the pro- 
file flattened, and converged with the gravity-only model 
near R V i r - After correcting for the gas fraction, the 
entropy profile was close to the self-similar predictions 
of gravitational structure formation simulations, as ob- 
served in massive galaxy clusters. 

5. In contrast with the Virgo and Perseus clusters, 
there was no evidence of further flattening in the entropy 
profile outside ^Rsooi nor unusually high f gas measure- 
ments. This indicates that significant gas clumping can- 
not be ubiquitous near ~R v i r , at least for relaxed (fossil) 
groups. 

6. There is good agreement between the best-fitting 
model for the gas emission and the Chandra, Suzaku and 
Rosat data, as well as the agreement in the cosmic X- 
ray background (the dominant background component) 
measurements made with each satellite. 

7. The dominant uncertainty in the gas properties at 
large radius, the shape CXB spectrum, can be signifi- 
cantly improved upon by using a Chandra mosaic of the 
outermost Suzaku regions to resolve the CXB into its 
component point sources so its properties can be con- 
strained directly. 
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